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Chapter  I 
INTRODUCTION 

Gas-liquid  contacting  with  a  continuous  liquid  phase  is  an  essential 
process  in  chemical  and  biochemical  industries.   Many  reactions  are 
controlled  by  the  transfer  of  a  component  from  the  gas  phase  to  the 
liquid  phase.   The  mass  transfer  characteristics  of  gas-liquid  contacting 
devices,  therefore,  have  been  extensively  studied. 

Gas-liquid  contacting  is  traditionally  provided  by  bubble  columns  and 
by  agitated  vessels.   Many  researchers  have  shown  that  separating  the  bubble 
column  into  stages  results  in  improved  mass  transfer  characteristics. 
Similar  results  have  been  presented  for  the  stage-wise  separation  of 
agitated  vessels. 

One  of  the  more  promising  contactor  designs  for  high  mass  transfer  rates 
with  low  power  consumption  is  the  airlift  tower.   The  airlift  tower  is  a 
bubble  column  divided  into  two  sections.   The  gas  is  sparged  into  one 
section,  causing  the  dispersion  density  in  that  section  to  be  lower  than 
the  density  in  the  unsparged  section.   The  resulting  pressure  difference 
causes  liquid  circulation.   An  upward  liquid  velocity  is  observed  in  the 
sparged  section,  and  a  downward  liquid  velocity  is  observed  in  the  unsparged 
section.   The  downward  liquid  flow  entrains  some  gas,  therefore  mass  transfer 
can  take  place  in  all  parts  of  the  tower. 

A  discussion  of  the  effects  of  material,  process  related  and  geometric 
parameters  on  mass  transfer  is  presented  in  Chapter  II.   The  mass  transfer 
coefficients  reported  for  the  airlift  tower,  bubble  column,  agitated  vessel, 
staged  bubble  column  and  staged  agitated  vessel  are  compared. 

Staging  has  been  shown  to  improve  the  mass  transfer  characteristics 
of  bubble  columns  and  agitated  vessels.   It  is  likely,  therefore,  the 


staged  airlift  tower  will  exhibit  superior  mass  transfer  characteristics 
as  compared  to  the  single  stage  airlift  tower.   A  comparison  of  one-stage 
and  two-stage  airlift  towers  is  presented  in  Chapter  III.   The  comparison 
is  based  on  mass  transfer  rates  and  efficiencies  at  superficial  gas 
velocities  ranging  from  0.120  m/s  to  0.455  m/s. 

The  optimal  design  and  operation  of  large  scale  airlift  towers  are 
contingent  upon  the  successful  modeling  and  prediction  of  performance. 
An  understanding  of  the  gas  phase  flow  is  necessary  to  model  the  airlift 
tower.   The  gas  phase  in  the  sparged  region  (upflow)  is  present  as  bubbles 
in  turbulent  upward  flow.   The  bubble  flow  in  the  unsparged  (downflow) 
region  is  more  complex.   The  forces  acting  on  a  bubble  in  the  downflow 
region  are  the  buoyant  force,  acting  in  the  upward  direction,  and  the 
drag  force,  acting  in  the  downward  direction.   The  bubble  flow  in  the 
downflow  region,  therefore,  has  both  upward  and  downward  components  of 
flow. 

In  Chapter  IV  a  statistical  method  is  introduced  to  analyze  the  bubble 
flow  for  an  air-water  system  in  the  downflow  section  of  a  one-stage  airlift 
tower.   This  method  is  theoretically  developed,  experimentally  verified  and 
applied  to  the  airlift  tower  to  measure  local  time  averaged  bubble  velocities 
and  void  fractions.   The  flowrates  of  the  gas  phase  in  the  upward  and  down- 
ward direction  were  determined  at  two  axial  positions  for  superficial  gas 
velocities  ranging  from  0.0682  m/s  to  0.3115  m/s.   The  bulk  liquid  circulation 
rate  and  regional  hold-ups  were  obtained  to  characterize  the  airlift  tower. 
Chapter  V  presents  conclusions  and  recommendations  for  future  work. 


Chapter  II 
OXYGEN  TRANSFER  IN  FERMENTATION  SYSTEMS 
Introduction 


Submerged  aerobic  fermentation  has  found  application  in  areas  ranging 
from  waste  treatment  to  manufacture  of  pharmaceuticals  and  food  products. 

Of  great  interest  is  the  production  of  single  cell  protein  (SCP)  for 

(1-3) 
animal  or  human  consumption     .   Substrates  for  SCP  production  include 

energy  source  materials  such  as  petroleum  derivatives,  waste  material 

such  as  bagasse  and  animal  manure,  and  renewable  source  material  such 

(4) 
as  cellulose  and  sugar    .   The  economic  viability  of  production  of  SCP 

from  these  substrates  hinges  in  large  degree  upon  the  growth  rate  of 

the  microorganisms  employed. 

Of  the  nutrients  necessary  for  microbial  growth,  oxygen  has  the  lowest 
solubility  in  the  liquid  phase.   The  design  of  fermentation  equipment  is 
oriented,  therefore,  toward  maximizing  the  oxygen  supply  available  for  cell 
growth  while  minimizing  the  power  requirement    . 

Schlegel    has  recently  reported  that  oxygen  can  be  successfully 
supplied  to  a  culture  by  adding  hydrogen  peroxide  to  a  shaking  flask 
under  nitrogen  blanket.   Cell  concentrations  achieved  were  comparable  to 
those  achieved  in  shaking  flasks  under  open  atmosphere.   Most  fermentations, 
however,  depend  on  the  mass  transfer  of  oxygen  from  the  gas  phase  to  the 
liquid  to  provide  the  required  oxygen. 

The  fermentation  process  typically  entails  a  three  phase  system, 
containing  the  cell,  medium  and  gas  phases.   The  mechanism  of  oxygen 
transfer  in  such  a  system  involves  three  steps  in  series.   The  oxygen  is 
transferred  from  the  gas  phase  to  the  liquid  phase,  and  then  from  the 
liquid  phase  to  the  cells.   At  the  cell-liquid  interface  a  chemical  reaction 
takes  place  with  the  oxidizable  substrates  through  the  respiratory  enzymes 


( 7  R*l 
of  the  cell.   '    The  oxygen  transfer  steps  may  encounter  a  gas-film 

resistance  and  a  liquid  film  resistance  at  the  gas-liquid  interface,  a 

liquid  film  resistance  at  the  liquid-cell  interface,  as  well  as  a  bulk 

(9) 
liquid  diffusional  resistance   . 

A  detailed  treatment  of  the  transport  equations  for  multi-phase  systems 

has  been  presented  by  Whitaker    .   The  treatment  associates  every  point 

in  space  with  an  averaging  volume  large  enough  to  give  rise  to  a  continuous, 

well-behaved,  space-averaged  value  of  a  state  variable  $  for  each  phase. 

For  momentum  transport,  9  represents  the  momentum  vector,  pv.   Whitaker's 

time  smoothed  relationship,  modified  by  Gray     and  applied  to  mass  transport 

in  a  given  phase  "i",  is 

3C. 

■5J-  +  v.  ■  Kt   =  V  •  (D.  •  VC.)  -  Z  K.  a.^C.  -  C.°)  +  R.     (1) 

where  C  is  the  molar  concentration,  v  is  the  mass  average  velocity,  t  is 
time,  D  is  the  mass  dispersion  tensor,  K  is  the  overall  mass  transfer 
coefficient,  a  is  the  interfacial  area  between  phases  i  and  j,  C°  is  the 
equilibrium  molar  concentration,  and  R  is  the  molar  rate  of  production 
due  to  reaction.   This  relationship  holds  under  the  assumptions  that  no 
phase  change  occurs,  the  fluid  is  incompressible,  the  concentration  of 
the  diffusing  species  is  small  and  the  fluctuation  terms  resulting  from 
time  averaging  are  negligible. 

A  fermentation  system  can  be  described  by  equation  (1)  with  oxygen  as 
the  diffusing  species.  When  the  oxygen  transfer  resistance  of  the  cell  is 
negligible,  the  cell  can  be  treated  as  a  reacting  species.  The  production 
term  in  equation  (.1)  is  replaced  by  the  volumetric  oxygen  consumption  rate 
of  the  cells  as  follows: 


ac. 

■ft-  +  v.  •  VC.  =  V  •  (D  •  VC.)  -  I   K.  a.j  (C±  -   C.°)  -  rX    (2) 

where  r  is  the  specific  oxygen  uptake  of  the  cells  and  X  is  the  cell  concen- 
tration.  If  the  fermentor  is  well  nixed  the  concentration  gradients  are 
zero.   Thus,  we  have 

vr  -  J  Ki  aij  (<V  -  ci)  - rx  ») 

The  gas  film  resistence  has  been  shown  to  be  negligible  in  fermentation 
systems       ,  and  thus,  in  the  liquid  phase  the  relationship  becomes 
3C 

it =  h*  (c* _  c)  ' rX  w 

where  K^a  is  the  overall  oxygen  transfer  coefficient  based  on  the  liquid  side 
resistance  controlling,  and  C*  is  the  phase  equilibrium  dissolved  oxygen  concentratioi 

Mass  transfer  in  a  gas-liquid  contacting  system  is  affected  by  a  number 
of  parameters  which  can  be  classified  as  material,  process  related  or 
geometric. 

Material  Parameters 

For  low  liquid  viscosity  systems  the  interfacial  properties  are  the 
important  material  parameters.   In  slow  bubble  formation  at  an  open  vertical 
tube  immersed  in  a  liquid  the  bubble  grows  until  its  buoyancy  exceeds  the 
surface  tension  holding  it  to  the  tube.   The  surface  tension  in  this  case 
controls  the  bubble  size  and  hence,  the  interfacial  area.   In  rapid  bubble 
formation,  however,  the  effect  of  surface  tension  is  less  pronounced,  and 
for  constant  flow  conditions  with  a  large  pressure  drop  across  the  sparger 


the  surface  tension  becomes  insignificant  compared  to  the  inertial  forces. 

In  pure  water,  bubbles  will  coalesce  easily  after  being  sparged,  but 
the  addition  of  impurities  hinders  coalescence.   This  has  been  explained 


(15) 


in  terms  of  a  change  in  surface  viscosity  or  elasticity  due  to  the  addition 
of  surface  active  compounds,   '  but  such  compounds  are  not  the  only  sub- 
stances which  affect  this  observation.   The  addition  of  electrolytes  to 
water  hinders  coalescence  because  of  the  formation  of  electrical  double  layers 
around  gas  bubbles  which  become  negatively  charged  and  repel  each  other. 
Zlokarnik     has  shown  that  a  salt  concentration  of  10  g/1  enhances  the 
absorption  rate  five-fold  and  its  influence  increases  until  at  100  g/1 
the  enhancement  is  almost  seven-fold.   The  surface  pressure,  defined  as 
the  difference  between  the  surface  tension  of  water  and  that  of  the  solution, 
was  shown  to  be  an  insignificant  parameter  in  the  coalescence  phenomena. 

A  surface  tension  gradient,  caused  by  a  temperature  or  concentration 
gradient  along  the  interface,  can  create  surface  turbulence  and  increase 
the  interfacial  area.   This,  in  turn,  increases  the  mass  transfer  rates    . 
Interfacial  properties  seem  to  affect  mass  transfer  primarily  through 
the  interfacial  area  by  promoting  or  hindering  coalescence  and  interfacial 
turbulence. 

Process  Related  Parameters 

The  mass  transfer  coefficients  in  gas-liquid  systems  have  been  empiri- 
cally or  semi-empirically  related  to  process  related  parameters.   The 
parameters  employed  include  the  volume  flowrate  of  gas  per  unit  liquid 
volume,  q/V,  and  the  power  consumption  per  unit  liquid  volume,  P/V. 

Frequently  the  power  consumed  in  mixing  under  ungassed  conditions,  P  , 

o 

is  employed. 

The  oxygen  transfer  coefficient  in  a  simple  bubble  column  has  been 
related  to  the  superficial  gas  velocity  by  an  equation  of  the  form 

V  =  l  vs  (5) 


where  C  and  a  are  constants  that  depend  on  column  geometry,  sparger  design 

.   ...    .   (18) 
and  void  fraction 

In  a  bubble  column  with  mechanical  agitation  the  correlation  is  made 
in  terms  of  the  power  of  agitation 

V  ■  *  C  (6) 
where  I  and  n  are  constants.   This  works  well  for  the  special  case  where 
the  power  consumed  in  aeration  is  small.   Miller     and  Tojo     have 
recommended  consideration  of  the  effective  power  as 

V  =   P  +  I  P  (7) 
e    m    2  a  v  ' 

where  P  is  the  effective  power  consumption,  P   is  the  power  consumed  in 
8  m 

mixing  and  P  is  the  power  consumed  in  compression  of  the  gas.   I„  is  a 
weighting  factor.   The  oxygen  transfer  coefficient  now  becomes 

V  ■  h    (Pm  +  h   Pa'"  (8) 
where  I  and  n  are  constants  that  depend  on  the  type  of  contacting  device. 
The  Miller  correlation  can  be  used  to  describe  systems  not  employing 
mechanical  agitation  by  assigning  a  zero  value  to  P  ,  making  it  a  general 
correlation. 

The  power  consumption  can  be  related  to  material  and  system  parameters 

such  as  the  Weber  number,  Aeration  number  and  liquid  and  dispersion  densities. 

Such  a  correlation  is  presented  by  Hassan  and  Robinson 

p  PT 

T0  -    *  V  NA   (")  (9) 

where  I  ,  and  n  are  functions  of  irapellor  size,  tank  size  and  electrolytic 
nature  of  the  aqueous  phase. 


Geometric  Parameters 

The  primary  geometric  consideration  is  the  type  of  contacting  device. 
Some  of  the  fermentors  are  shown  in  Tables  1  and  2.   They  are  classified 
as  single  and  multi-stage  fermentors. 

Single  Stage  Fermentors 

The  fermentors  used  traditionally  in  industry  are  the  bubble  column  and 
agitated  tank,  generally  operated  in  a  batch  mode.   The  chief  advantage  of 
bubble  columns  is  that  they  have  no  moving  parts,  thus  eliminating  the  need 
for  seals  and  reducing  maintenance  expense.   Bubble  columns  require  a  small 
amount  of  floor  space,  and  fermentation  heat  transfer  requirements  are 
easily  provided  through  the  tower  walls.   Mashelkar(22)  has  reported  that 
coalescence  will  cause  a  reduction  in  interfacial  area  which  reduces  the 
oxygen  transfer  rate  in  bubble  columns  with  length-to-diameter  ratios  greater 
than  12.   This  effect  is  dependent  upon  the  material  parameters. 

Yoshida  et  al.     have  reported  an  oxygen  transfer  coefficient  of 

6.7  min   at  a  superficial  gas  velocity  of  1333  cm/min  for  a  sodium  sulfite- 

(23) 
air  system.   Malshelkar     has  found  a  mass  transfer  coefficient  of  7.8  min"1 

at  a  superficial  gas  velocity  of  2400  cm/min  for  absorption  of  CO   in  a 

sodium  carbonate,  sodium  bicarbonate  buffer  solution.   Other  representative 

values  of  mass  transfer  rates (   '25)  are  shown  in  Table  1. 

The  efficiency  of  oxygen  transfer  can  be  shown  by  the  performance 

ratio  defined  as  the  rate  of  oxygen  transfer  divided  by  the  power  consumed 

in  achieving  that  transfer  rate.    ' 

C  =  P  (10) 

A  high  value  of  the  performance  ratio  indicates  a  high  oxygen  transfer 
efficiency.   The  bubble  column  typically  exhibits  a  high  performance  ratio 


at  low  oxygen  transfer  rates,  but  as  shown  in  Fig.  1,  this  value  decreases 
rapidly  with  an  increase  in  the  oxygen  transfer  rate.   This  is  due  primarily 
to  bubble  coalescence  and  inadequate  liquid  phase  mixing. 

To  increase  the  oxygen  transfer  rates,  agitation  is  applied  to  the 
bubble  column.   This  agitation  increases  liquid  mixing  and  the  interfacial 
area.   The  liquid  and  gas  phases  can  be  considered  as  completely  mixed, 
and  since  a  high  air  flow  rate  is  not  needed  to  generate  mixing,  the  gas 
input  to  the  agitated  fermentor  can  be  reduced  considerably. 

Fukuda  et  al.     reports  oxygen  transfer  rates  ranging  from  2.6  to  13.8  min 

in  a  turbine  agitated  sodium  sulfite-air  system.  Their  results  are  reviewed  by 

(27) 
Hospodka  et  al.      Agitated  tanks  are  characterized  by  higher  oxygen 

transfer  rates  than  found  in  bubble  columns,  but  the  total  power  consumption 

is  high.   The  performance  ratio  can  be  seen  in  Fig.  1  to  be  low  for  agitated 

tanks . 

Many  novel  gas-liquid  contacting  devices  have  been  investigated, 

including  several  that  depend  on  gas  entrainment  to  achieve  oxygen 

/TO   o  y_  TTN 

transfer   '      .   One  of  the  more  promising  of  these  designs  is  the 
air-lift  tower  invented  by  Lefrancois  et  al.     in  1954.   This  device, 
shown  in  Fig.  2,  can  be  described  as  a  bubble  column  divided  into  two 
sections.   Air  is  sparged  into  one  of  these  sections,  causing  the  disper- 
sion in  that  section  to  have  a  lower  density  than  the  ungassed  section. 
The  resulting  pressure  difference  causes  liquid  circulation.   In  essence, 
the  airlift  tower  provides  liquid  agitation  and  mixing  by  using  a  high 
gas  flowrate.   The  airflow  rate  in  the  airlift  tower  is  generally  larger 
than  that  in  the  agitated  tank,  and  the  need  for  mechanical  agitation  is 
eliminated.   The  airlift  tower  has  been  described  by  Hatch     as  having 


three  regions:   the  upflow  region,  the  downflow  region,  and  the  head  region 
consisting  of  the  two  phase  dispersion  above  the  upflow  and  downflow  regions. 
It  is  in  the  head  region  that  gas  disengagement  takes  place.   Gas  bubbles 
are  entrained  by  the  downward  flow  of  liquid  in  the  downflow  region  so 
oxygen  transfer  can  take  place  in  all  parts  of  the  airlift  tower. 

The  airlift  tower  has  no  moving  parts;  therefore,  the  need  for  seals 
is  eliminated.   Maintenance  expenses  are  low  as  compared  to  the  agitated 
tank,  and  the  downward  flow  of  liquid  from  the  head  region  of  the  tower 
reduces  some  of  the  foaming  problems  found  in  fermentation  in  bubble  columns 
at  similar  oxygen  transfer  rates.   The  performance  of  the  airlift  tower  in 
mass  transfer  applications  is  good. 

(19) 
Hatch     found  oxygen  transfer  rates  of  up  to  6.3  min"1  for  a  concentric 

cylinder  airlift  with  Candida  intermedia  on  n-alkanes.   Gasner(32)  investigated 

a  rectangular  airlift  with  a  baffle  separating  the  upflow  and  downflow  sections 

and  found  a  l^a  of  9.7  min"1.   He  has  suggested  that  the  rectangular  system 

inherently  has  less  coalescence  than  the  concentric  cylinder  airlift  due  to 

configurational  differences.   A  further  modification  of  the  airlift  tower 

is  described  by  Belfield(33)  who  found  a  l^a  of  7.8  min"1  in  a  cylinder 

split  into  two  equal  sections  by  a  baffle.   This  split  cylinder  airlift 

has  a  lower  wetted-surface-to-volume  ratio  than  the  concentric  cylinder 

airlift  tower  which  increases  the  liquid  circulation  rate  and  hence  the 

liquid  mixing.   The  ease  of  construction  of  the  split  cylinder  airlift  makes 

this  fermentor  design  attractive. 

(35) 
Ho     has  modeled  the  airlift  tower  by  treating  the  upflow  region  as 

N  tanks- in-series,  the  downflow  region  as  M  tanks-in-series  and  the  head 

region  as  one  completely  mixed  tank  in  which  gas  disengagement  takes  place. 
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This  approach  can  be  modified  to  incorporate  a  dispersion  model.  The  use 
of  such  detailed  models,  however,  depends  upon  further  characterization 
of  the  airlift  tower.   The  amount  of  gas  entrained  in  the  downflow  region 
as  a  function  of  incoming  air  flowrate,  axial  position  and  tower  design 
should  be  established. 

Multi-Stage  Fermentors 

A  number  of  motivations  exist  for  considering  multi-stage  operation  in 
a  fermentation  system.   Chen  et  al.(36)  and  Moo-Young(37)  have  pointed  out 
the  value  of  partial  mixing  in  continuous  fermentation.   Chen  et  al/37' 
found  a  single  stirred  tank  followed  by  a  plug  flow  reactor  to  be  optimal 
in  an  autocatalytic  fermentation  with  sterile  feed  and  a  high  desired 
conversion  of  substrate  to  cell  mass.   Several  authors(38_42)  have  found 
multi-stage  systems  useful  in  continuous  fermentations  where  more  than 
one  biological  reaction  take  place.   In  batch  and  continuous  operations 
improved  oxygen  transfer  rates  have  frequently  been  found  in  multi-stage 
operation. 

Various  researchers    '   ■  5>43-45)  have  found  gQod  oxyggn  transfer 
characteristics  in  bubble  towers  filled  with  ceramic  or  screen  packing.   A 
mass  transfer  coefficient  of  7.8  min"1  has  been  found  by  Voyer  and  Miller(45) 
for  a  C02-H20  system  in  a  tower  packed  with  open  end  cylindrical  mesh 
packing  (0.5  inch  x  0.5  inch)  at  a  superficial  gas  velocity  of  3658  cm/min. 
The  use  of  perforated  plates  in  multistage  sieve  tray  bubble  towers  has 
been  shown  to  improve  mass  transfer  characteristics  as  compared  to  the 
single  stage  bubble  column(25'45'47) .   Hsu  et  al.(25'46)  have  compared  the 
multistage  sieve  tray  tower  to  a  tower  packed  with  Koch  motionless  mixers 
and  found  the  Koch  mixers  to  improve  the  oxygen  transfer  characteristics. 
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Oxygen  transfer  coefficients  of  3.33  min~  were  found  at  a  superficial  gas 
velocity  of  180  cm/min  for  a  mixed  culture  fermentation. 

Falch  and  Gaden   '    and  Paca  and  Gregr(50'51)  have  added  mechanical 
agitation  to  each  stage  of  the  sieve  tray  bubble  tower  to  further  improve 
the  oxygen  transfer  rates.   This  type  of  fermentor  could  be  considered  as 
a  multistage  agitated  tank  tower.   Paca  and  Gregr^50^  have  found  an  oxygen 
transfer  coefficient  of  18.6  min"  at  a  superficial  gas  velocity  of  89.3 
cm/min  and  an  agitation  speed  of  800  rpm.   Another  multi-stage  tower 
investigated  was  the  vibrating  disk  tower  by  Tojo  et  al/52^  for  low  oxygen 
transfer  rates. 

Calderbank   '    and  Mashelkar   }    have  shown  that  coalescence  will 
occur  in  tall  bubble  columns,  which  in  turn  reduces  the  interfacial  area. 
Orazem     has  shown  that  an  optimal  height  exists  for  oxygen  transfer  in 
a  split  cylinder  airlift  tower  using  the  sulfite  oxidation  method.   This 
may  also  be  a  result  of  bubble  coalescence.   The  concept  of  employing  a 
multi-stage  airlift  tower  to  enhance  oxygen  transfer  appears  to  be  valid. 


13 

Notation 

a 

interfacial  area 

A 

cross-sectional  area 

C 

Molar  concentration 

C° 

Equilibrium  molar  concentration 

C* 

Saturation  liquid  phase  oxygen  concentration 

d 

impellor  diameter 

D 

Mass  dispersion  tensor 

I 

Constant 

K 

Overall  mass  transfer  coefficient 

V 

Liquid  side  oxygen  transfer  coefficient 

m 

constant 

n 

constant 

N 

impellor  rotational  speed 

NA 

3 
Aeration  Number,  Q/Nd 

X 

N 
We 

Oxygen  transfer  rate 

2  3 
Weber  Number,  N  d  p  la 

P 

Total  power  consumption 

P 
a 

Power  consumed  in  aeration 

Pe 

Effective  power  consumption 

P 

IE 

Power  consumed  in  mixing  in  a  gassed  system 

po 

Power  consumed  in  mixing  in  an  ungassed  system 

Q 

Volumetric  gas  sparge  rate 

r 

Specific  oxygen  uptake  of  the  cells 

R 

molar  rate  of  production  due  to  reaction 

t 

Time 

V 

Velocity 
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v  Superficial  gas  velocity,  Q/A 

X       Cell  concentration 


Subscripts 

L  Liquid 

D  Dispersion 

Greek 

a  Constant 

K  Performance  ratio 

p  Density 

a  Surface  tension 
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.      Performance  Ratio  as  a  Function  of  Oxygen  Transfer  Rate; 
1.      rectangualr  airlift    (32),    2.      one  stage  split   cylinder 
airlift    (33),    3.      concentric  cylinder  airlift    (18),    4.   agitated 
tank   (26),    5.      aeration  tower    (5). 
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OXYGEN  TRANSFER  RATES  AND  EFFICIENCIES  IN  ONE  AND  TWO  STAGE  AIRLIFT  TOWERS 

Introduction 


Submerged  aerobic  fermentation  has  found  application  in  areas  ranging 
from  waste  treatment  to  manufacture  of  pharmaceuticals  and  food  products. 

Of  great  interest  is  the  production  of  single  cell  protein  (SCP)  for  animal 

fX— 31 

or  human  consumption      .   Substrates  considered  for  SCP  production 

include  energy  source  materials  such  as  petroleum  derivatives,  waste 

materials  such  as  bagasse  and  animal  manure,  and  renewable  source  material 

(4) 
such  as  cellulose  and  sugar    .   The  economic  viability  of  production  of 

SCP  from  these  substrates  hinges  in  large  degree  upon  the  growth  rate  of 

the  microorganisms  employed. 

Of  the  nutrients  necessary  for  the  growth  of  microorganisms,  oxygen 

has  the  lowest  solubility  in  the  liquid  phase.   The  design  of  fermentation 

equipment  is  geared,  therefore,  toward  maximizing  the  oxygen  supply  available 


for  cell  growth  while  minimizing  the  power  requirement 


(5) 


Schlegel    recently  reported  that  oxygen  can  be  successfully  supplied 
to  a  culture  by  adding  hydrogen  peroxide  to  a  shaking  flask 
under  nitrogen  blanket.  Cell  concentrations  achieved  were  comparable  to 
those  achieved  in  shaking  flasks  under  open  atmosphere.  Most  fermenta- 
tion equipment,  however,  depend  on  the  mass  transfer  of  oxygen  from  the> 
gas  phase  to  the  liquid.   The  bubble  column  and  the  agitated  tank  are 
traditional  gas-liquid  contacting  devices. 

High  oxygen  transfer  rates  were  found  by  Hatch'  '   '    in  a  concentric 
cylinder  airlift  tower  with  low  power  requirements.   The  airlift  tower 

provides  liquid  agitation  through  the  use  of  a  large  gas  flow  rate. 

(8) 

Belfield    described  an  airlift  tower  which  consists  of  a  cylinder 

divided  into  two  sections  by  a  baffle.   The  reduced  wetted  surface  area  to 
volume  ratio  and  its  ease  of  construction  make  the  split  cylinder  airlift 
tower  attractive. 
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Multistage  systems  have  been  shown  to  have  good  mass  transfer  char- 

(9-13) 
acteristics.   Various  researchers       have  found  improved  mass  transfer 

characteristics  in  bubble  towers  filled  with  random  ceramic  or  screen 

packing.   Mass  transfer  coefficients  of  7.80  min   were  found  by  Voyer 

(13) 
and  Miller     for  a  CO  -  H.O  system  in  a  tower  packed  with  open  end 

cylindrical  mesh  packing  (0.5  inch  x  0.5  inch)  at  a  superficial  gas  velocity 

of  3658  cm/min.   The  use  of  perforated  plates  in  multistage  sieve  tray 

bubble  towers  was  shown  to  improve  mass  transfer  characteristics  as  compared 

to  the  single  stage  bubble  column       .   Hsu  et  al.    '    compared  the 

multistage  sieve  tray  tower  to  a  tower  packed  with  Koch  motionless  mixers 

and  found  the  Koch  mixers  to  improve  the  oxygen  transfer  characteristics. 

Oxygen  transfer  coefficients  of  3.33  min   were  found  at  a  superficial  gas 

velocity  of  180  cm/min  for  a  mixed  culture  fermentation. 

Falch  and  Gaden   '    and  Paca  and  Gregr     added  mechanical  agitation 
to  each  stage  of  the  sieve  tray  bubble  tower  to  further  improve  the  oxygen 
transfer  rates.   This  type  of  fermentor  could  be  considered  as  a  multistage 
agitated  tank  tower.   Paca  et  al. (19)  found  an  oxygen  transfer  coefficient  of 
18.6  min   at  a  superficial  gas  velocity  of  89.3  cm/min  and  an  agitation 
speed  of  800  rpm  for  a  sulfite-air  system.   Another  multistage  tower  investi- 
gated was  the  vibrating  disk  tower  by  Tojo  et  al. 

In  an  investigation  preceding  the  work  presented  here,  the  effect  of 

column  height  on  oxygen  transfer  in  a  split  cylinder  airlift  tower  was 

(21) 
studied  for  a  sodium  sulfite-air  system    .   An  optimal  baffle  height 

for  a  15.24  cm  diameter  column  was  found  to  be  near  60  cm  for  a  baffle  placed 

6.20  cm  above  the  base.   It  was  concluded  that  a  multi-stage  airlift  tower 

would  merit  further  consideration  as  a  means  of  increasing  oxygen  transfer 

rates. 
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The  objective  of  this  study  was  to  determine  the  applicability 
of  a  multistage  split  cylinder  airlift  tower  for  gas-liquid  mass  transfer 
such  as  expected  in  fermentation  systems.  The  range  of  superficial  gas  velo- 
cities studied  was  735  cra/min  to  2728  cm/min.   The  volume  of  gas  per  volume 
liquid  per  minute  (wm)  ranged  from  3.28  min   to  12.21  min  .  This  range 
was  chosen  to  provide  good  data  for  scale-up.   The  volumetric  gas  flow 
rate  can  be  described  as  the  product  of  the  superficial  gas  velocity  and 
the  cross-sectional  area. 

Q  =  vs  •  A  (1) 

The  wm  ratio  is  seen  to  be 

wm  -  2_  (2) 

where  h  is  the  tower  height .   Thus 

v 
wm  -  —  (3) 

h 

Industrial  fermentations  are  typically  conducted  with  a  wm  ratio  of 
approximately  three.   In  order  to  keep  superficial  velocities  constant 
in  tower  scale-up  the  studies  at  low  tower  height  must  be  conducted  at 
high  wm  values.   Thus  studies  in  a  one-meter  tall  tower  with  a  wm  of 
twelve  can  be  related  to  a  four-meter  tall  tower  with  a  wm  ratio  of  three. 

Materials  and  Methods 

Two  split  cylinder  airlift  towers  were  employed  in  this  study.   They 
are  shown  in  Fig.  1.   Each  tower  held  a  liquid  volume  of  20.5  I   and  had 
an  inside  diameter  of  15.24  cm.   The  single  stage  tower  had  a  baffle  height 
of  107.1  cm  and  the  gap  between  the  lower  edge  of  the  baffle  and  the  bottom 
of  the  column  was  5.5  cm.   The  cross  sectional  area  allowed  for  flow  between 
the  bottom  of  the  column  and  the  baffle  was  equal  to  the  upflow  and  downflow 
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cross  sectional  areas.   Air  was  sparged  into  the  upflow  section  of  the 
column  through  a  0.35  cm  I.D.  plexiglas  tube  which  was  level  with  the 
lower  edge  of  the  baffle. 

The  two  stage  airlift  was  composed  of  two  sections,  each  with  a  50.8  cm 
high  baffle.   The  distance  between  the  lower  edge  of  the  baffle  and  the 
bottom  of  the  tower  section  was  5.5  cm,  as  in  the  single  stage  tower.   The 
two  sections  were  divided  by  a  plate  with  a  sparger  identical  to  the  one 
used  in  the  single  stage  tower.   The  sections  were  designed  so  as  to  allow 
a  large  range  of  superficial  gas  velocities  without  restricting  the  head 
region  volume  for  the  lower  stage.   To  provide  a  level  control  for  the  upper 
stage  an  external  downcomer  was  placed  at  an  adjustable  height  at  the  top 
of  the  head  region.   The  lower  stage  level  control  was  provided  by  pumping 
continuously  from  the  desired  liquid  level  to  the  upper  stage.   The  total 
circulation  between  the  two  stages  provided  in  this  manner  was  27  ml/min 
for  a  volume  exchange  time  of  6.3  hours.   Air  was  filtered  and  the  flowrate 
was  monitered  by  a  Fisher  and  Porter  (T2-1308/2)  rotameter  which  had  been 
calibrated  by  use  of  a  wet  test  meter. 

Open  end  manometers  were  used  to  obtain  pressure  drop  and  local  hold 
up  data.   These  manometers  were  located  at  the  bottom  and  the  top  of  each 
baffle  on  the  upflow  and  downflow  sections  of  the  towers.   The  overall 
gassed  liquid  level  was  measured  to  determine  the  overall  gas  holdup. 

The  oxygen  transfer  characteristics  were  found  through  use  of  the 
sodium  sulfite  oxidation  method  by  Cooper  et  al.("^.   This  method  can 
be  safely  employed  in  a  direct  comparison  of  fermentation  equipment.   A  one 
weight  percent  solution  of  technical  grade  sodium  sulfite  was  used  with  one 

part  per  million  cobaltous  catalyst.   The  measured  solution  pH  was  10.0 

(23) 
as  recommended  by  Fari  and  Wang 
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The  sodium  sulfite  and  the  catalyst  CoSO,  were  allowed  to 
dissolve  thoroughly  in  a  large  tank.   The  solution  was  then 
pumped  into  the  tower  being  used.   The  sparging  of  air  at  the  desired 
flowrate  began  at  a  recorded  time.   Fifty  milliliter  samples  were 

withdrawn  at  regular  intervals  and  from  each  sample  a  10  ml  aliquot  was 

(24) 
titrated  using  the  iodometric  analysis  for  sodium  sulfite  concentration 

The  reaction  of  sodium  sulfite  to  form  sodium  sulfate  is  zero  order  with 
respect  to  sodium  sulfite  under  the  conditions  used,  so  the  zero  order 
reaction  rate  constant  could  be  used  to  determine  the  amount  of  oxygen 
transferred  under  pseudo  steady  state  conditions  with  regard  to  the  liquid 
oxygen  concentration.   The  liquid  oxygen  concentration  was  assumed  to  be 
negligible  and  the  apparent  ICa  was  calculated  assuming  perfect  liquid 
mixing . 

The  comparative  oxygen  transfer  characteristics  of  the  one  stage  and 
two  stage  towers  were  found  for  gas  superficial  velocities  ranging  from 
720  cm/min  to  2728  cm/min.   The  liquid  was  at  an  average  temperature  of  18°C 

and  the  literature  values  used  for  the  saturated  oxygen  concentration  were 

,    c  (25) 

adjusted  for  temperature  variation 
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Results  and  Discussion 

Gas  Hold-up 

The  overall  gas  hold-up  was  determined  for  the  two-stage  tower,  the  single 
stage  tower,  and  the  single  stage  tower  packed  with  three  Koch  LY  static 
mixers  and  packed  with  seven  Koch  LY  static  mixers.   The  static  mixers 
consisted  of  three  one-inch  corrugated  plastic  sheets  glued  together  and 
cut  so  as  to  fit  inside  the  upflow  section  of  the  airlift  tower.   The 
direction  of  mixing  was  tangential  to  the  baffle.   The  Koch  mixers     are 
described  by  Hsu  et  al.     and  Hsu  et  al.       The  overall  gas  hold-up 
is  shown  in  Fig.  2.   The  values  for  the  two-stage  tower  are  the  arithmetic 
average  values  for  both  stages.   The  gas  hold-up  in  the  single  stage  tower 
is  seen  to  be  much  lower  than  in  the  two  stage  tower  for  high  superficial 
gas  velocities.   At  a  superficial  gas  velocity  of  2262  cm/min  the  gas  hold-up 
in  the  two  stage  tower  is  40  percent  higher  than  in  the  one  stage  tower. 
This  could  be  attributed  to  increased  coalescence  in  the  single  stage  tower, 
causing  the  bubbles  to  increase  in  size  and  escape  more  quickly  from  the 
dispersion.   A  reduced  gas  residence  time  would  result. 

The  gas  holdup  in  the  head  region,  upflow  region  and  downflow  region 
are  shown  in  Fig.  3  for  the  upper  stage  of  the  two-stage  tower  and  in 
Fig.  A  for  the  single  stage  tower.   Results  for  the  lower  stage  are 
similar  to  those  shown  in  Fig.  3.   The  very  high  hold-up  in  the  head  region 

of  the  two  stage  tower,  0.575  at  a  superficial  gas  velocity  of  2778  cm/min, 

(281 
could  be  partly  due  to  entrainment  of  surface  air.   Calderbank  .   .  found 

surface  aeration  to  be  important  in  agitated  gas  liquid  contacting  systems. 

The  entrainment  of  surface  air  due  to  the  intense  turbulence,  as 

well  as  the  downward  velocity  of  the  liquid,  could  play  an  important  role 

in  mass  transfer. 
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The  hold-up  in  the  head  region  of  the  single  stage  tower  was  not  as 
high  as  that  in  the  two  stage  tower.   At  a  superficial  gas  velocity  of 
2262  cm/min  the  observed  hold-up  was  0.355.   The  agitation  in  the  single  stage 
tower  was  observed  not  to  be  as  intense  as  in  the  two  stage  tower,  and  the 
carry-through  of  entrained  air  in  the  downflow  section  was  not  as  large. 
Coalescence  and  the  resulting  rapid  escape  from  the  system  would  account 
for  these  observations. 

Oxygen  Transfer  Characteristics 

Under  the  assumptions  that  the  content  of  the  aqueous  phase  is  completely 
mixed  and  the  oxygen  content  in  gaseous  phase  remains  invariant  with  respect 
to  the  spatial  position,  an  oxygen  balance  over  the  aqueous  phase  in  each 
stage  of  the  airlift  tower  gives 

£  .  -Na  +  V  (C*  -  C)  (4) 

Under   a  pseudo-steady  state  condition,    this   reduces   to 

NA  =  Kja    (C*  -  C)  (5) 

The  La  to  be  determined  here  is  on  a  liquid  volume  basis,  i.e.,  the  concen- 
trations are  in  units  of  milligrams  0  per  unit  liquid  volume 
and  not  per  unit  dispersion  volume. 

NA 

V  -  cTTT  (« 

The  saturated  liquid  oxygen  concentration  value  was  found  from  the  literature 
for  a  1  wt  percent  salt  solution  at  the  dispersion  temperature  and  atmospheric 
pressure.   The  effect  of  pressure  on  the  saturated  liquid  oxygen  concentration 
was  neglected  for  the  short  towers  used,  since  the  pressure  recorded  at  the 
bottom  of  .the  towers  was  close  to  atmospheric.   The  liquid  oxygen  con- 
centration was  assumed  to  be  zero.   This  assumption  of  a ~ maximum  oxygen 
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transfer  driving  force  leads  to  a  conservative  estimate  of  the  apparent 
oxygen  transfer  coefficient. 

The  oxygen  transfer  coefficient,  ICa,  was  based  on  the  liquid  film 

resistance  controlling.  This  assumption  is  reasonable  for  oxygen  transfer 

(28,29) 

in  gas-liquid  contacting  systems 

ILa  is  shown  in  Fig.  5  as  a  function  of  superficial  gas  velocity 
for  the  two-stage  and  one-stage  towers.   At  a  superficial  gas  velocity  of 
2262  cm/min  the  IC  a  was  54  percent  higher  in  the  two-stage  tower  than  in 
the  single  stage.   Reproducibility  of  the  data  was  within  a  standard  devia- 
tion of  six  percent.   The  rise  in  La  begins  to  fall  off  at  high  gas  velocities 
for  the  one  stage  tower.   This  is  probably  due  to  the  effect  of 
coalescence,  and,  while  IC  a  is  a  linear  function  of  superficial  gas 
velocity  for  the  two  stage  tower  at  the  velocities  studied,  a  leveling 
can  be  expected  at  some  yet  undetermined  higher  gas  velocity.   The  tL  a 
found  in  the  lower  stage  in  the  two-stage  tower  matched  consistently  that 
found  in  the  upper  stage  within  experimental  accuracy.   This  could  be 
due  to  the  relatively  small  consumption  of  the  available  oxygen  in  the 
sparged  air.   The  calculated  oxygen  mole  fraction  in  the  exit  gas  stream 
at  a  superficial  gas  velocity  of  2400  cm/min  was  0.198.   That  in  the  gas 
entering  the  second  stage  0.204.   An  inlet  oxygen  mole  fraction  of 
0.21  was  assumed  in  these  calculations. 

The  oxygen  transfer  coefficient  is  shown  in  Fig.  6  for  the  single 
stage  tower  packed  with  three  Koch  LY  static  mixers  and  with  seven  Koch  LY 
static  mixers.   The  effect  of  the  mixing  elements  was  to  reduce  the  liquid 
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circulation  rate  and  to  decrease  the  oxygen  transfer.   Under  the  conditions 
investigated,  the  empty  airlift  tower  was  found  to  have  better  oxygen 
transfer  characteristics  than  the  airlift  tower  with  packing  in  the  upflow 
section.   The  static  mixers  did  increase  the  interfacial  area  by  reducing 
coalescence  but  this  effect  was  overshadowed  by  the  decrease  in  liquid 
velocities. 

It  must  be  noted  that  the  static  mixers  employed  were  not  of  optimal 
configuration.   The  mixing  was  not  three-dimensional  and  the  corrugation 
was  coarse  relative  to  the  column  diameter.   These  results  should  not  be 
interpreted  to  imply  that  packing,  or  Koch  mixers,  are  not  desirable 
under  any  conditions  in  air-lift  or  split  cylinder  air-lift  towers. 

As  samples  were  withdrawn  from  the  tower  during  the  course  of  a  run, 
the  total  liquid  volume  decreased  by  two  percent,  corresponding  to  an 
ungassed  liquid  level  change  of  2.2  cm  in  the  one-stage  tower.   The  effect 
on  K_ a  of  varying  the  ungassed  liquid  level  from  9  cm  above  the 
baffle  to  9  cm  below  is  shown  in  Fig.  7.   The  effect  of  a  variation  of 
2  cm  in  liquid  level  can  be  seen  to  be  negligible. 

Oxygen  transfer  is  enhanced  in  the  two  stage  tower  as  compared  to  the 
single  stage  at  large  superf icialgas  velocities.   An  important  consideration 
in  the  evaluation  of  the  data  presented  is  the  total  volume  occupied  by 
the  aeration  equipment.   The  oxygen  transfer  characteristics  were  adjusted 
to  a  dispersion  volume  basis  by  the  relation 

Kja'  =  i^a  (1-H)  (7) 

where  t^a'  is  the  oxygen  transfer  coefficient  on  a  dispersion  volume  basis 
and  H  is  the  overall  gas  hold-up.  The  comparison  of  oxygen  transfer  coeffi- 
cients on  a  dispersion  volume  basis  is  more  reasonable  than  a  similar 
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comparison  on  a  tower  volume  basis.  The  two-stage  tower  was  designed  to 
allow  a  large  void  volume  between  the  first  and  second  stages.   This  void 
volume  can  be  greatly  reduced  in  tower  design  for  a  specific  gas  flowrate, 
while  the  volume  occupied  by  the  liquid-gas  dispersion  cannot  be  reduced 
without  changing  the  system  flow  characteristics. 

The  oxygen  transfer  coefficient  based  on  a  dispersion  volume  basis 
is  shown  in  Figure  8  for  the  two-stage  and  one-stage  towers.   The  ICa' 
values  for  the  two  stage  tower  still  exceed  those  of  the  one-stage  tower  at 
high  superficial  gas  velocities  but  the  ICa'  is  not  a  linear  function  of  gas 
velocity  as  it  was  for  the  liquid  volume  based  ICa.   At  2000  cm/min  the 
K^a'  rises  less  rapidly  with  gas  velocity.   The  ICa'  based  on  dispersion 
volume  was  27  percent  higher  in  the  two-stage  tower  than  in  the  one  stage 
tower  at  a  superficial  gas  velocity  of  2262  cm/min. 
Performance  Ratio 

As  stated  in  the  introductory  section,  the  design  of  fermentation 

equipment  hinges  on  maximizing  the  oxygen  transfer  rates  while  minimizing 

the  power  requirements.   The  performance  ratio  defined  by  Hatch    is  a 

measure  of  the  oxygen  transferred  perunit  of  power  required  to  transfer  that 

oxygen.   The  adiabatic  power  of  compression  for  air  is 

P  0.286 
P  =  0.0153  Q^ffcr)      -  1]  (8) 

1 

where  P  is  the  power  requirement,  hp,  Q  is  the  gas  flowrate,  SCFM,  P  is 

the  atmospheric  pressure,  and  ?2   is  the  pressure  at  the  gas  sparger,  psia. 

The  power  term  obtained  is  characteristic  of  the  tower  only  and  not  of 

the  air  line  between  the  compressor  and  the  tower. 
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The  performance  ratio,    ?,    is 

s  =  yp  do) 

The  performance  ratio  is  shown  in  Fig.  9  for  the  one-  and 
two-stage  towers.   The  towers  are  comparable  at  low  oxygen 
transfer  rates  but  at  a  transfer  rate  of  180  m  mole  0./Z  hr  the  two-stage 
tower  becomes  clearly  more  efficient.   The  performance  ratio  has  its  highest 
values  at  low  oxygen  transfer  rates  in  the  two-stage  tower  and  an  optimal 
value  near  150  m  mole  0,/jl  hr  for  the  one-stage  tower. 

Figure  10  displays  the  performance  ratio  for  the  packed  one-stage 
airlift  tower.   The  results  for  the  two  stage  tower  are  also  shown  in  this 
figure  for  reference. 

Comparison  with  Literature 

The  oxygen  transfer  coefficient  as  a  function  of  superficial  gas  velocity 
is  presented  in  Fig.  11  for  various  airlift  tower  systems.   Values  were 
adjusted  to  a  liquid  temperature  of  18°C.   The  physical  characteristics 
and  dimensions  of  the  various  towers  studied  are  shown  in  Table  1.   It 
is  important  to  note  that  while  our  work  and  those  of  Belfield    and  Gasner   ' 
utilized  a  sulfite-air  system,  the  work  presented  by  Hatch  ,   made  use  of  a 
fermentation  system.   Thus,  the  IL  a  values  obtained  by  Hatch  are  not 
directly  comparable  to  those  obtained  through  the  use  of  the  sulfite 
oxidation  method. 

Both  curves  2  and  4  represent  single  stage  split  cylinder  airlift 
towers.   Flooding,  or  the  appearance  of  a  coalescence  effect  occurs  in 
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Belfield's  tower  at  a  superficial  gas  velocity  of  600  cm/ sec  while  in  the 
smaller  tower  it  occurs  at  roughly  1200  cm/sec.   This  discrepancy  may  be 
due  in  part  to  the  lack  of  an  optimal  baffle  height  in  the  larger  tower. 

The  results  from  Gasner  clearly  show  the  merit  of  using  a  wide  baffle 
in  an  airlift  system.   Large  oxygen  transfer  coefficients  were  obtained 
at  low  gas  velocities.   This  may  have  been  due  in  part  to  the  use  of  a 
multiple  hole  sparger  as  compared  to  the  single  open  tube  sparger  used  in 
Belfield's  and  our  work.   The  use  of  a  similar  sparger  or  a  perforated 
plate  should  be  considered  in  the  split  cylinder  airlift.   The  rectangular 
airlift  should  be  investigated  at  higher  gas  velocities  and  in  a  multiple 
stage  combination. 

The  performance  ratios  are  shown  in  Fig.  12  for  various,  systems.   The 
airlift  towers  are  clearly  superior  to  the  bubble  column  and  agitated  tank 
fermentors.   The  application  of  the  multi-stage  concept  to  the  split 
cylinder  airlift  tower  yields  a  high  efficiency  at  high  oxygen  transfer 
rates.   The  performance  ratio  of  a  small  system  will  be  smaller  than 
that  of  a  large  system  due  to  a  larger  surface-area-to-volume  ratio  of  the 
small  system.   Therefore,  the  performance  ratios  of  lines  4  and  5  may  be 
expected  to  increase  if  scaled  up  to  the  size  of  systems  1,  2  or  3. 
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Conclusions 

The  use  of  a  multi-stage  airlift  tower  has  merit  in  gas-liquid  contacting. 
Oxygen  transfer  coefficients  (£L  a)  of  over  25  min   (or  1500  hr   )  were  ob- 
tained at  a  superficial  gas  velocity  of  2728  cm/min  for  a  sodium  sulfite-air 
system  in  a  two-stage  tower.   At  a  superficial  gas  velocity  of  2262  cm/min 
the  YL  a  for  the  two-stage  tower  was  54  percent  higher  than  in  the  one-stage 
system. 

The  use  of  a  multi-stage  tower  requires  more  reactor  volume  than  the 
single  stage  tower  of  equal  liquid  volume,  but  the  oxygen  transfer  coeffi- 
cient based  on  the  total  dispersion  volume  is  still  27  percent  higher  than 
in  the  one  stage  tower  at  a  superficial  gas  velocity  of  2262  cm/min. 

At  low  gas  velocities,  the  two-stage  and  one-stage  towers  have  equiva- 
lent mass  transfer  characteristics.   At  1300  cm/min  the  1C  a  begins  to  level 
off  for  the  one-stage  system  while  the  fC  a  for  the  two  stage  tower  increases 
linearly  with  superficial  gas  velocity  in  the  range  studied. 

The  two-stage  tower  becomes  considerably  more  efficient  than  the  single 
stage  tower  at  oxygen  transfer  rates  larger  than  180  m  mole/8,  hr. 
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Notation 

2 
A       cross  sectional  area,  cm 

C  liquid  oxygen  concentration,  mg/£ 

H  gas  holdup  based  on  dispersion  volume 

h  tower  height,  cm 

5L  a  oxygen  transfer  coefficient  based  on  liquid  volume,  min 

ICa'  oxygen  transfer  coefficient  based  on  dispersion  volume,  min 

N  rate  of  oxygen  consumption,  mg/£  min 

P  adiabatic  power  of  compression,  hp 

Q  volumetric  gas  flowrate,  SCFM 

t  time,  min 

Vq  superficial  gas  velocity,  cm/rain 

wm  volume  of  gas  per  volume  of  liquid  per  minute,  min 

Superscript 

*       saturation 

Greek  letter 

£       performance  ratio,  lb  0„/hp  hr   (or)   g  0„/J 
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LIQUID    VOLUME  20.5.      f 

DIAMETER  15.2    Cm 

SPARGER  I.D.        0.95  Cm 


Fig.    1.      Schematic  Diagram  of   Single   Stage  and   Two   Stage  Tower. 
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Fig.  2.   Overall  Gas  Hold-Up  as  a  Function  of  Superficial  Gas 
Velocity;  O  Two  Stage  Tower, A  One  Stage  Tower,  O  One 
Stage  Tower  with  Threo  Koch  LY  Mixers, G  One  Stage 
Tower  with  7  Koch  LY  Mixers. 
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3.   Gas  Holdup  as  a  Function  of  Superficial  Gas  Velocity 
for  a  Two-Stage  Airlift  Tower  in  the  Upper  Stage; 
&   Head  Region,  O  Upf low  Region,  C  Downflow  Region. 
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APPENDIX  A:   Sulfite  Oxidation  Technique 
The  use  of  sulfite  oxidation  to  characterize  the  oxygen  transfer  of 
aeration  equipment  was  proposed  by  Cooper  et  al.     in  1944.   This  method 
utilizes  the  very  rapid  reaction  of  sulfite  with  oxygen  in  the  presence 
of  copper  or  cobalt  catalyst  to  form  sulfate.   The  oxygen  consumption  may 
be  determined  directly  through  analysis  of  incoming  and  exiting  gas 
streams  for  oxygen  or  through  reaction  stochiometry.   Although  the 

reaction  has  been  studied  extensively  for  the  last  two  decades,  consid- 

(2-5) 
erable  confusion  still  exists  as  to  the  reaction  kinetics     .   This 

confusion  seems  to  be  due  to  the  apparent  complexity  of  the  reaction 

kinetics.   According  to  Astarita    et  al. ,  when  the  reaction  is 

catalyzed  by  cobalt,  it  is  zero  order  with  respect  to  oxygen  at  a  sulfite 

concentration  of  0.06  M,  first  order  at  0.25  M,  and  second  order  at  0.25  M 

to  1.0  M.   The  reaction  is  zero  order  with  respect  to  sulfite  for  a  5  wt% 

(2) 
Na2S03  solution^   . 

(2) 
Fan  and  Wang    have  found  that  the  reaction  proceeds  most  rapidly  in 

a  pH  range  of  9  to  11,  and  the  reaction  rate  is  fastest  at  a  temperature  of 

90°C.   The  existence  of  an  optimal  temperature  was  attributed  to  the 

increase  of  reaction  rate  with  temperature  coupled  with  a  decrease  of 

oxygen  solubility  in  the  sulfite  solution  with  increasing  temperature. 

They  would  not  identify  explicitly  the  dependence  of  the  reaction  rate 


on  the  cobalt  catalyst  concentration  from  1  ppm  to  10  ppm  CoSO  .   Fan  and 

(2)  -! 

Wang    report  reaction  rate  constants  exceeding  0.063  gmole  Na,S0,/(m  s) , 

(7)  , 

and  Orazem  and  Erickson    report  reaction  rates  up  to  0.182  gmole  Na.SO  /(m  s) 

using  1  ppm  cobaltous  catalyst.   Robinson  and  Engel    and  Finn    have 

recommended  the  use  of  cobalt  as  a  catalyst  over  copper.   Pirt  et  al. 

found  that  oxygen  absorption  rates  were  two  to  five  times  higher  when  the 
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reaction  was  catalyzed  by  cobaltous  ions  than  when  it  was  catalyzed  by 
Cupric  ions. 

Maxon  and  Johnson     claimed  that  the  absorption  of  oxygen  by  a 

sulfite  solution  catalyzed  by  copper  must  be  gas  film  controlled.   Schultz 

(12) 
and  Gaden  "  studied  oxygen  transfer  through  a  stagnant  air-sulfite 

solution  interface.   Agitation  of  the  gas  phase  did  not  increase  the 
oxygen  transfer  rate;  therefore,  they  concluded  that  the  gas-film  was 
not  controlling. 

Under  the  assumptions  that  the  content  of  the  liquid  phase  is  com- 
pletely mixed  and  the  oxygen  content  of  the  gaseous  phase  remains  invariant 
with  respect  to  the  spatial  position,  an  oxygen  balance  over  the  aqueous 
phase  yields 

§  -  -NA  +  KLa  (C*  -  C)  (1) 

where  NA  is  the  rate  of  oxygen  consumption  by  reaction,  K  a  is  the  oxygen 
transfer  coefficient,  C*  is  the  saturation  liquid  phase  oxygen  concentration 
and  C  is  the  liquid  phase  oxygen  concentration.   Since  the  oxidation  of 
sulfite  is  very  rapid,  the  reaction  is  mass  transfer  controlled,  and  the 
bulk  oxygen  concentration  can  be  assumed  negligible. 

Thus,  a  pseudo  steady  state  condition  is  attained,  in  which  the  rate 
of  oxygen  consumption  by  reaction  is  balanced  by  the  rate  of  oxygen  transfer 
from  the  gas  phase  to  the  liquid,  i.e., 

NA  "  V  C*  (2) 

The  oxygen  transfer  coefficient,  l^a,  can  be  found  by  measuring  the  rate 

of  oxygen  consumption. 

NA 
V  =  C^  (3) 
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The  sulfite  oxidation  reaction  is  zero  order  with  respect  to  the  sodium 
sulfite  concentration.   Thus,  equation  (3)  becomes 

V  _  _  c*  fit  w 

where  AC   s()  is  the  change  in  sodium  sulfite  concentration,  mg  Ha  SO./t, 
in  the  time  interval  At,  min,  C*  is  in  units  mg  0.,/i,  and  K  is  the  stoichio- 
metric conversion  factor,  0.127  mg  Cymg/mg  Na^Cy   The  value  of  C*  can 
be  found  in  the  literature,  or  can  be  determined  with  a  dissolved  oxygen 
probe.   The  rate  of  oxygen  consumption  can  also  be  determined  through 
comparison  of  the  mole  fraction  of  oxygen  in  the  inlet  and  exit  gas  streams. 

The  value  of  t^a  determined  through  sulfite  oxidation  usually  differs 
from  values  found  under  fermentation  conditions.   Greenhalgh  et  al.(13) 
found  that  for  low  cell  concentrations  the  oxygen  transfer  coefficient 
under  fermentation  conditions  was  lower  than  the  oxygen  transfer  coefficient 
determined  through  sulfite  oxidation.   At  high  cell  concentrations  the 
l^a  under  fermentation  conditions  was  higher  than  that  determined  through 
sulfite  oxidation.   The  presence  of  microbial  cells  as  a  solid  phase  has 
been  shown  to  enhance  oxygen  transfer       . 

Benedek  and  Bennett     note  that  the  f^a  found  through  sulfite  oxida- 
tion is  not  equivalent  to  that  found  for  an  air-water  system.   This  is 
attributed  to  the  difference  in  liquid  properties  between  the  sulfite 
solution  and  water  and  to  the  enhancement  of  mass  transfer  by  the  reaction. 
The  addition  of  electrolytes  to  water  results  in  the  formation  of  electrical 
double  layers  around  gas  bubbles  which  become  negatively  charged  and  repel 
each  other.   This  effect  prevents  coalescence  and  can  also  significantly 
increase  oxygen  transfer.   Zlokarnik(18)  reports  a  seven-fold  increase  in 
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ILa  for  a  5  wt.  percent  Na  SO,  solution  when  compared  to  water.   The 
sulfite  method,  therefore,  can  be  used  safely  only  as  a  means  of  comparing 
the  relative  oxygen  transfer  efficiencies  of  two  systems. 

In  utilizing  the  sulfite  oxidation  technique  for  a  comparison  of 
oxygen  transfer  devices,  a  number  of  precautions  should  be  followed. 

1.  The  solution  properties  must  be  constant  from  experiment  to 
experiment.   The  addition  of  a  sulfite  salt,  as  noted  by  Benedek  and 
Bennett    ,  can  greatly  change  the  viscosity,  density,  surface  charge  and 
interfacial  tension  of  the  liquid-gas  dispersion.   A  change  in  salt  concen- 
tration will  also  affect  the  saturation  value  of  the  dissolved  oxygen  con- 
centration, which  then  changes  the  driving  force  for  mass  transfer. 

If  not  stored  in  an  airtight  container,  solid  sulfite  will  slowly  react 
with  the  oxygen  in  the  air  to  form  sulfate.   Constant  initial  sodium  sulfite 
concentrations  in  repeated  experiments  can  be  attained  by  adding  increasing 
amounts  of  solids  to  compensate  for  the  sulfite  degradation.   This,  however, 
makes  comparison  from  experiment  to  experiment  difficult  since  the  fluid 
and  interfacial  properties  are  highly  dependent  upon  salt  concentration   '    . 
Since  the  zero  order  reaction  rate  constant  is  to  be  determined,  the  main- 
tenance of  a  constant  initial  sulfite  concentration  is  not  necessary.   The 
total  salt  concentration  is  the  parameter  which  should  be  held  constant  for 
all  experiments  to  be  carried  out  for  comparisons. 

2.  The  assumption  of  a  negligible  liquid  oxygen  concentration  results 
in  a  maximum  oxygen  transfer  driving  force  and  the  iC  a  determined  will  be 
conservative.   If  an  actual  dissolved  oxygen  concentration  exists,  the 

measured  ICa  will  be  low.   It  is  important  to  allow  the  reaction  to  proceed 

(2) 
under  near  optimal  conditions.   The  pH  should  be  between  9  and  11  '  and  a 


57 


cobaltous  catalyst  should  be  employed.   Liquid  mixing  also  should  be  such 
that  no  strong  sulfite  concentration  gradient  exists  as  this  could  result, 
in  a  non-zero  oxygen  concentration  at  some  point  in  the  contactor. 

The  assumption  of  a  negligible  oxygen  concentration  in  the  liquid 
phase  can  be  checked  by  the  use  of  a  dissolved  oxygen  probe. 

3.   Liquid  samples  should  be  analyzed  as  quickly  as  possible  after 
being  withdrawn  from  the  contactor.   Samples  should  be  stabilized  immedi- 
ately, as  described  by  "Standard  Methods  for  the  Examination  of  Water  and 
Wastewater" (19). 
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Table  Bl.  General  Tower  Characteristics 


UP  upflow  section 

DN  downflow  section 

L   lower  edge  of  baffle  (5.5  cm   above  bottom  of  stage) 

U  upper  edge  of  baffle 


1 

One  St 

age 

Two  Stage 

Liquid  Volume 

20.  5C 

20.50 

Diameter 

15.2C 

cm 

15.20 

cm 

Baffle  Height 

107. 0C 

cm 

(2)  51.00 

cm 

Distance  from  Stage  Bottom 

to  Bottom  Edge  of  Baffle 

5.5C 

cm 

(2)  5.50 

cm 

Sparger  Tube  Height 

5.5C 

cm 

(2)  5.50 

cm 

Sparger  Tube  I.  D. 

0.9! 

cm9 

(2)  0.95 

cm, 

Upflow  Cross-sectional  Area 

1 

90. 7C 

cm" 

90.70 

cm 

Pressure  Tap  Number 

Baffle 

Baffle 
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and  Location 
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Edge 

Section 

Edge 
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Table  B5.  Pressure  Data  for  the  One-Stage  Airlift  Tower  -  Unpacked 


Run 
No. 


9 
10 
11 

25 
26 
27 


Tap  No. 


Pressure  (cm  HO  gauge) 


0 


104.3 
104.3 
104.8 
106.4 
105.3 
104.  S 
104.8 
105.1 


114.5 
120.1 
127.7 
123.4 
122.6 
127.7 
128.7 
130.2 


106.4 
106.1 
106.9 
107.4 
106.4 
106.4 
106.4 
106.9 


112.5 
117.5 
124.6 
122.1 
121.8 
125.2 
126.7 
127.7 
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Height 

(cm) 


30.5 
33.5 
36.2 
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[Air  Pressure  (psig) 

Air  Flow  Rate  (std.  m3/min) 

Superficial  Gas  Velocity  (cm/min) 

vvm  (min"'-) 

Liquid  Temperature  (°C) 

Salt  Concentration  (wt  %) 

Reaction  Rate  (mg  Na?S0-,/£  min) 

KLa  (min-1)            J 

Power  Consumption  (W) 

Oxygen  Transfer  Rate  (g  02/«  hr) 

Performance  Ratio  (g  02/J  x  10^) 

Hold  Up: 
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Head 
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Chapter  4 
STATISTICAL  ANALYSIS  OF  BUBBLE  FLOW  IN  THE  DOWNFLOW  SECTION  OF  THE  AIRLIFT  TOWER 

Introduction 

Gas-liquid  contacting  with  a  continuous  liquid  phase  is  an  important 
operation  in  chemical  and  biochemical  industries.   The  bubble  column  and  stirred 
tank  are  traditionally  employed  in  gas-liquid  mass  transfer.   One  of  the 
more  promising  contactor  designs  for  high  mass-transfer  applications  is  the 
airlift  tower  invented  by  Lefrancois  et  al. (1)   in  1954.   This  device,  shown 
in  Fig.  1,  can  be  described  as  a  bubble  column  divided  into  two  sections. 
Air  is  sparged  into  one  of  these  sections,  causing  the  dispersion  in  that 
section  to  have  a  lower  density  than  that  in  the  ungassed  section.  The 
resulting  pressure  difference  causes  a  liquid  circulation.  Liquid  agitation 
is  provided,  therefore,  by  the  gas  supply. 

The  airlift   tower  has  been  described  by  Hatch    as  having  three 
regions:   the  upflow,  downflow  and  head  regions.   The  head  region, 
where  gas  disengagement  takes  place,  is  the  two  phase  dispersion 
above  the  upflow  and  downflow  regions.   Gas  bubbles  are  entrained  by  the 
liquid  flow  in  the  downflow  region;  therefore,  both  phases  are  present  in  all 
parts  of  the  airlift  tower. 

(3) 
Ho  et  al.     has  modeled  the  airlift  tower  by  treating  the  upflow  region  as 

N  tanks-in-series,  the  downflow  region  as  M  tanks-in-series  and  the  head 

region  as  one  completely  mixed  tank.   This  approach  can  be  modified  to 

incorporate  dispersion  models.  The  gas  phase  in  the  upflow  region  consists 

of  bubbles  in  turbulent  upward  flow.   The  flow  in  the  downflow  section  of 

the  tower  is  more  complex.   The  forces  on  a  bubble  in  the  downflow  section 

are  the  buoyant  force  and  the  drag  force,  acting  in  opposite  directions. 

The  net  direction  of  flow  for  a  large  bubble  will  be  upward  and  that  for  a 

small  bubble  will  be  downward.   There  are,  therefore,  significant  upward 
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and  downward  components  of  gas  flow  in  the  downflow  section  of  the  airlift  tower. 
The  proportion  of  the  incoming  gas  entrained  in  the  downflow  section  of  the 
airlift  tower  needs  to  be  established  as  a  function  of  spatial  position, 

incoming  gas  flow  rate  and  system  parameters  for  use  in  detailed  models,  such 

(3) 
as  that  proposed  by  Ho  et  al.    . 

The  two-phase  flow  in  the  airlift  tower  is  turbulent,  and,  therefore, 

a  stochastic  analysis  is  often  more  appropriate  than  a  deterministic  approach 

(4) 
for  the  airlift  tower    .   Analyses  of  turbulence  have  long  used  statistical 

techniques,  employing  time  smoothing  to  obtain  mean  values  of  state  variables, 
autocorrelation  to  identify  the  cyclical  nature  of  a  process,  and  cross- 
correlation  to  find  the  time  lag  between  dependent  signals  generated  by  the 
turbulent  flow   '   . 

Mesch  et  al.    describe  the  use  of  correlation  methods  to  measure  the 
velocity  of  moving  surfaces.   Ong  and  Beck    found  bulk  slurry  velocities 
in  turbulent  slurry  flow  by  finding  the  cross  correlation  between  signals  from 
pressure  transducers  placed  a  known  distance  apart.   The  signals  were  generated 
by  the  pressure  variation  inherent  in  turbulent  flow.   Ong  and  Beck    re- 
ported a  measurement  accuracy  of  two  percent.   Beck  et  al.    and  Mesch 
and  Kipphan     measured  bulk  particle  velocities  in  pneumatic  conveyors 
by  obtaining  the  cross  correlation  of  naturally  occurring  flow  noise,  measured 
by  two  capacitance  transducers  placed  a  known  distance  apart  along  the  axis  of  the 
conveyor.   Mesch  and  Kipphan     also  employed  optical  techniques.   Lassahn 
measured  the  flow  velocities  of  two-phase  fluids  by  using  radiation  beams  as 
sensors  to  detect  the  passing  of  inhomogeneities.   Fourier  spectrum  analysis 
of  the  resulting  signals  was  used  to  extract  the  velocity  information. 
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Cross-correlation  techniques  may  be  used  to  extract  time  lag  information 
from  two  related  signals.   The  method  used  to  generate  these  signals  determines 
the  interpretation  of  the  results  as  a  local  or  bulk  phenomenon.  Oki  et  al.(12) 
developed  an  optical  probe  which  could  be  inserted  into  a  bed  to  find  local 
solid  particle  velocities.  Herringe  and  Davis (13)  used  the  cross-correlation 
technique  to  determine  local  bubble  velocities  in  air-water  mixtures.  The 
signals  were  generated  by  two  resistivity  probes.  A  single  resistivity 
probe  will  yield  local  void  fractions  and  bubble  sizes  as  functions  of 
time       .  The  time  averaged  value  is  used  to  characterize  these  parameters. 

Cross-correlation  techniques  have  been  shown  to  be  useful  for  velocity 
measurement  in  stationary  stochastic  systems  with  a  single  average  velocity. 
Haines  et  al.     have  presented  bimodal  cross-correlation  functions,  indicating 
the  presence  of  two  average  solid  velocities  in  a  fluidized  bed.   The  cross- 
correlation  of  signals  obtained  from  the  downflow  section  of  the  airlift 
tower  can  be  expected,  therefore,  to  be  bimodal  since  two  significant  velocity 
components  are  present. 

The  intent  of  this  work  was  to  develop  a  technique  capable  of  extracting 
the  time  lag  and  the  time  fraction  of  flow  associated  with  each  velocity 
component  from  a  bimodal  cross-correlation  function.  A  further  objective 
of  the  work  was  to  experimentally  verify  the  validity  of  the  technique,  and 
to  demonstrate  its  use  in  the  complex  two-phase  flow  of  the  downflow  section 
of  the  airlift  tower. 

Theoretical 

A  point  in  a  two-phase  flow  field  can  only  be  occupied  by  one  phase 
at  a  time.  At  any  time,  t,  therefore,  the  void  fraction  at  a  point  in  a 
gas-liquid  system  can  only  have  a  value  of  unity  or  zero.   The  void  fraction 
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at  a  point  as  a  function  of  time  will  take  the  form  of  a  square  wave 
with  values  of  one  or  zero    .   The  expectation,  or  mean,  of  this 
function  is  the  time-averaged  void  fraction. 

Each  point  in  the  gas-liquid  flow  field  has  associated  with  it  a 
void  fraction  as  a  function  of  time.   Suppose  that  two  points  are  parallel 
to  the  net  direction  of  flow  and  are  separated  by  a  short  distance  so 
that  a  bubble  passing  through  one  point  can  reasonably  be  expected  to 
pass  through  the  other.   Then,  the  respective  void  fractions  as  functions 
of  time  will  be  related,  and  the  cross  correlation  between  these  two 
functions  will  give  rise  to  the  time  lag  between  them.   If  the  distance 
between  the  two  points  is  known,  the  bubble  velocity  can  be  determined. 

Suppose  that  the  gas  phase  flow  in  the  gas-liquid  flow  field  has 
significant  upflow  and  downflow  components.   If  point  X  is  above  point  Y, 
a  bubble  moving  downward  will  pass  first  through  X  and  then  Y.   Similarly, 
a  rising  bubble  will  pass  first  through  Y  and  then  through  X.  The  cross 
correlation  between  the  void  fraction  at  X  and  that  at  Y  is  bimodal.   Each 
peak  of  the  cross-correlation  function  is  characterized  by  a  time  lag. 
The  time  lag  corresponding  to  a  downward  velocity  is  positive  and  that 
corresponding  to  an  upward  velocity  is  negative.   A  time  fraction,  x, 
of  the  bubbles  passing  through  the  two  points  is  associated  with  the 
downward  velocity,  and  a  time  fraction  1-X,  is  associated  with  the  upward 
velocity. 

The  void  fraction  at  point  X  can  be  represented  as  the  time  series, 
X(t),  with  the  upward  component,  X  (t),and  the  downward  component 
Xdn(t),  i.e., 

xct)  =  xup(t)  +  xdn(t)  (1) 
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The  relationship  among  X(t),  X  (t)  and  X,  (t)  can  be  seen  in  Fie.  2. 

up        dn  ° 

Similarly,  the  void  fraction  at  Y  can  be  expressed  as 

Y(t)  =  Yup(t)  +  Ydn(t)  (2) 

When  a  bubble  is  present  at  point  X,  X(t)  has  a  value  of  one,  and  when 
the  liquid  phase  is  present,  X(t)  has  a  value  of  zero.   A  bubble  rising 
and  a  bubble  moving  downward  cannot  simultaneously  occupy  a  given  point; 
therefore,  the  upward  and  downward  components  cannot  simultaneously 
have  a  value  of  one.   The  components  can,  however,  simultaneously  have  a 
value  of  zero. 

The  sample  means  of  X(t)  and  Y(t)  converge  to  constant  values,  and 
therefore,  the  series  X(t)  and  Y(t)  can  be  considered  stationary.   The 
presence  of  bubbles  at  point  X  or  Y  can  be  considered  to  be  a  random 
occurrence.  The  time  series  representing  bubbles  moving  upward  at  X  is 
assumed  independent  of  the  time  series  representing  bubbles  moving  down- 
ward at  Y,  and  similarly,  the  series  representing  bubbles  moving  downward 
at  X  is  assumed  independent  of  the  series  representing  bubbles  moving 
upward  at  Y.   The  time  series  representing  the  bubbles  in  upward  flow  at 
points  X  and  Y,  however,  are  not  independent  since  the  bubbles  in  upward 
flow  will  pass  through  point  Y  before  passing  through  point  X.  These 
series,  X  (t)  and  Y   (t),  are  related  by  a  time  lag  which  can  be 

determined  through  correlation  techniques.   The  series  X,  (t)  and  Y   (t) 

dn         dn 

are  similarly  related  by  a  time  lag. 

The  cross-covariance  function,  Y   (x) ,  for  two  stationary  time 
series,  W(t)  and  Z(t),  is  a  function  of  time  lag,  x(17),  i.e., 

Ywz(x)  -  E{[W(t)  -  E[W(t)]]   [Z(t  +  i)  -  E[Z(t  +  x)]]}  (3) 
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Equation 

(3)  can 

be  expanded 

to  yield 

Y  (t) 
WZV  ; 

=  E{W(t)Z(t 

+  T)  - 

W(t)E[Z(t  +  t)]  -  Z(t  + 

r)E[W(t)] 

+  E[W(t 

)]  E[Z( 

t  +  T)]} 

(4) 

The  expectation 

of  a  sum  is 

equal 

to  the  sum  of  the  expectations    ; 

therefore, 

WT) 

=  E[W(t)Z(t 

+  T)] 

-  E[W(t)  E[Z(t  +  t)] 

-  E[Z(t 

+  t)  E[W(t)]]  +  E[E[W(t)]  E[Z(t 

+t)]]   (5) 

Since 

E[E[W(l 

:)]]   =  E[W(t 

01, 

(6) 

equation 

(5)  becomes 

V*> 

=  E[W(t)Z(t 

+  x)]  ■ 

-  E[W(t)]  E[Z(t  +  t)] 

(7) 

The 

cross-correlation  function 

is  found  by  normalizing  the  cross- 

covariance  function  as 

Pwz^ 

VT> 

(8) 

tV°>  V^ 

where  ¥,_ , 
WW 

(0)  and 

Yzz(0)  are, 

respectively,  the  autocovariance 

functions 

of  W(t)  and  Z(t) 

at  t  •  0.   The  autocovariance  function  of  W( 

t)  is 

given  by 

Yww^ 

=  E[W(t)W(t 

+  T)]  - 

■  E[W(t)]E[W(t  +  t)] 

(9) 

77 
and  the  autocorrelation  function  is  defined  as 

The  time  series  W(t)  can  be  discretized  to  form  the  discrete  series 
W(k)  where  k  is  an  integer  value  corresponding  to  a  sample  number.  Time,  t, 
is  the  product  of  the  sample  increment,  At,  and  the  sample  number,  k. 
Time  lag,  x,  can  be  expressed  as  the  product  of  the  sample  increment 
At  and  the  lag,  h,  defined  by 

h  --S- 

At 

The  cross-correlation  and  autocorrelation  functions  for  discrete  stationary 
series  are  functions  of  h.  An  estimator  of  the  autocorrelation  function 
of  W(k) ,  p  (h)  ,  is  given  by 

K-h 
Z   [(W(k)  -  W)  (W(k  +  h)  -  W)] 

Loo  =  ^ — ■ ■ (ID 


ww         K 


E   (W(k)  -  W) 
k=l 


where  K  is  the  total  number  of  observations,  and  W  is  the  estimated 
mean  of  W(t) .   The  estimator  of  the  autocorrelation  function  P   (h)  , 
approaches  the  true  autocorrelation  function  value  as  the  sample 
number  increases  according  to 

p   (h)  =  p   (h)  +  0  (K_l5).  (12) 

WW       WW 

The  calculation  of  the  autocorrelation  function  value  can  be  viewed 
as  the  matching  of  a  discrete  function  W(k)  with  another  function 
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identical  to  W(k) ,  but  shifted  by  some  lag,  h.   This  second  function  is 
W(k  +  h) .   Suppose  that  W(k)  is  a  discrete  series  of  random  pulses  with 
only  one  observation  for  every  pulse.   Clearly,  the  matching  at  a  zero 
lag  of  W(k)  with  W(k  +  0)  is  one  to  one.   As  shown  in  Fig.  3a,  the  match- 
ing is  not  one  to  one  for  W(k)  and  W(k  +  1) .  The  autocorrelation  function 
for  this  series  is  shown  to  have  a  value  of  one  at  h  =  0,  and  zero  at 
h  ?*  0  (Fig.  3b).   Suppose  that  the  function,  W(k)  ,  is  again  a  discrete 
series  of  random  pulses,  but  now  with  three  observations  for  every 
pulse.  Even  in  this  case  the  matching  of  W(k)  with  W(k  +  0)  is  one  to 
one;  however,  the  correlation  of  W(k)  with  W(k  +  1)  has  a  value,  as  does 
the  correlation  of  W(k)  with  W(k  +2),  as  shown  in  Fig.  4a.   The  auto- 
correlation of  W(k)  with  three  observations  for  every  pulse  still  has 
a  value  of  one  at  h  =  0,  but,  as  shown  in  Fig.  4b,  the  peak  is  spread. 

The  form  of  the  autocorrelation  and  the  cross-correlation  functions 
is  dependent  upon  the  number  of  observations  for  each  pulse,  I.      The 
autocorrelation  function  of  W(k)  is  intuitively  assumed  to  be  described  by 

£-1   'd-hi 
pWW(h)  =  (-~Y~)  <13a) 

where  d  is  the  theoretical  lag.   For  the  autocorrelation  function  the 
theoretical  lag  is  zero.   The  cross-correlation  function  is  assumed 
to  follow  equation  (13a) ,  but  in  this  case  the  theoretical  lag  may  be 
nonzero . 

t  -  1  |d_h| 
pwz(h)  =  (~T~)  C"b) 

The  empirical  justification  of  equations  (13a)  and  (13b)  will  be 
examined  later. 
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The  time  series  defined  in  equations  (1)  and  (2)  can  be  discretized 
to  become  functions  of  the  sample  number.   The  cross-covariance  function  is 
given  by  a  discrete  version  of  equation  (7)  which  is 


Vh)  ■  E[(Vk)  +  Xdn(k))  (Yup(k  +  «  +  Ydn(k  +  h))] 


E[Xup(W  +  Xdn(k)]  E[Yup(k  +  h)  +  Ydn(k  +  h)]    (14) 


Equation  (14)  can  be  expanded  and  rearranged  to  obtain 


V10  =  E[Xup(W  Yup(k  +  h)J  -  EtXup(k)]  E[Yup(k  +  h)] 


+  E[X  (k)Y,  (k  +  h)]  -  E[X   (k)]  E[Y,  (k  +  h)] 
up    an  up        an 


+  E[X.  (k)Y   (k  +  h)]  -  E[X   (k)]  E[Y   (k  +  h)  ] 
an    up  an        up 


+  E[Xdn(k)Ydn(k  +  h)]  -  E[Xdn(k)   E[Ydn(k  +  h)]   (15) 


Y   (h)  =  Y      (W  +  Y      (h)  +  Y      (h)  +  Y      (h)   (16) 

:  XY       X  Y         X  Y ,        X,  Y         X,  Y . 
up  up        up  dn        dn  up        dn  dn 

The  cross  covariance,  therefore,  is  additive.  Since  X   and  Y,  are 

up      dn 

assumed  independent,  as  are  X,   and  Y   ,  this  relationship  becomes 

dn      up 


YXY(h)  =  TX  Y   (h)  +  YX  Y   (h)  (17) 

up  up       dn  dn 
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From  equation  (  8 ) ,  the  cross-covariance  function  becomes 

Y   (h)  -  p       (h)[y       (0)  Y       W]h 

XY       X  Y       XX       Y  Y 

up  up      up  up      up  up 

+  P      (ti)  [Y      (0)  Y      W]h 

X,  Y,       X,  X,      Y,  Y,               (18) 
dn  dn       dn  dn      dn  dn 

By  the  nature  of  the  one-zero  square  wave  series,  we  have 

[X(k)]2  -  X(k)                                        (19) 

Thus,  from  equation  (9), 

Y      (0)  =*  U          -  U   2   ,                            (20) 

XX        XX 
up  up        up     up 

Y      (0)  -  V          -  u   2   ,                            (21) 

dn  dn        dn     dn 

Y      (0)  =  v          -   V        2   ,                            (22) 

Y  Y         Y      Y 
up  up       up     up 

Y      (0)  •  y    -  y   2                              (23) 

Y  Y        Y      Y 
dn  dn       dn     dn 

Introducing  equation  (13b)  and  equations  (20)  through  (23)  into  equation 

(18) 

the  cross-covariance  function  becomes 

,.  .    ,1   -  ljdup"h!  .,           2, ....           l.Ji 

y     (h)  -  (    i            [Cu   -  v        ).(u   -  v        )]2 

XY                    X      X       Y      Y 

up     up      up     up 

,  ,t  -   ljd,  -  h|  ,,                         2,  .           2,, h 
+  (  l     )   1  dn     1    [(ux   -  ux   )  (uy   -  uy   )] 

dn    '  dn      dn     dn 

(24) 
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where  d 

up 

represents 

the  theoretical  lag  corresponding  to  the  upward 

velocity  i 

component, 

and  d   is  the  theoretical  lag  corresponding  to  the 

downward  velocity  component.   The  lags,  d   and  d   ,  are  opposite  in 

sign  and, 

in  this  study,  have  a  magnitude  between  50  and  150  in  units 

of  h.  The  value  of 

the  cross-covariance  function  at  a  lag  of  d   is, 

up 

from  equation  (24) , 

f      (dup)  = 

Ly   "  y     )  (y     ~  y     )  ] 

XY  up 

XX       Y      Y 
up   up      up     up 

,  ,1  -   lj  dn"  up'  ,,                        2W         l^h   ,„,  , 
+  ( — jj— )         t(wx  -  ux   )(yy   uy   )]   (25a) 

dn    dn     dn   dn 

A  similar 

expression 

can  be  found  for  d.  . 
dn 

w  ■ 

,1   -  lj  up    dn  „           2S            2,  ,5j 

(  I     >'          KlJX   "  UX   )(UY   "  UY   )] 
up     up     up     up 

+  [(ux  -  yx  2)(uy  -  my  2)]h                           (25b) 

dn    dn     dn    dn 

T-he  cross- 

-covariance 

function  for  the  combined  upward  and  downward 

component 

will  have 

local  maxima  at  d   and  d,  ,  respectively,  just  as  the 
up      dn 

individual  cross-covariance  functions  for  unidirectional  upward  and  down- 

ward flow 

will  exhibit  maxima  at  d   and  d,  ,  respectively, 
up      dn 

Even 

though  the 

cross  covariance  at  h  =  0  is ,  from  equation  (24), 

Y„(0)  ■  (" 

I   -  1,   up  '  , ,           2W           2^  ,ij 
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it  is  reasonable  to  assign  a  value  of  zero  to  y   (0) .   The  reason  is 
that  zero  time  lag  corresponds  to  an  infinite  bubble  velocity,  which  is 
improbable  in  the  present  system.   The  cross-covariance  function> 

therefore,  is  linearly  translated  so  that  the  value  of  y      (0)  is  zero. 

The  value  of  the  cross-covariance  function  at  d   relative  to  f   (0) 

up  XY 

is 

L   =  Y   (d   )  -  y   (0)  ,  (27) 

up   TXY  up;    YXY^  '' 
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Similarly  at  d,  ,  we  have 
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The  peak   fraction  at   d        is   defined   as 
up 


L 


up        L        +  L , 
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(31) 
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Since  d   ,  d,   and  (d   -  A.   )  are  large  in  this  study,  we  have 
up   dn       up    dn 

f-   1)     '"up  "   "dn'    .   Q 

These  assumptions  will  be   experimentally  justified  later. 

From  equations    (28),    (30),    (31)    and    (32),    the  peak  fraction  becomes 

[(Ux  -      Ujj  Xily  -      Uy  )] 
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If  every  bubble  which  passes  through  point  X  also  passes  through  point  Y, 

u    =  u    =  m  (34) 

X     HT      up 
up     up 


and  similarly, 

\     =  *I.     '   Wdn  (35) 
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The  peak  fraction  at  d   can  then  be  expressed  as 
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When  u,      and  u        are  small,   equation    (37)    becomes 
an  up 

11 

Q       =  H2 


up       y       +  u 

up         an 


A  similar  expression  results  for  the  peak  fraction  at  d,   as 

dn 


(38a) 


dn   u   +  u 

on    up 
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As  mentioned  previously,  the  means,  u,   and  u   ,  are  respectively  the 

dn      up 

time  averaged  void  fractions  resulting  from  the  downward  and  upward 

components  of  flow.   Hence,  equations  (38a)  and  (38b)  become  respectively 

* 
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where  the  total  void  fraction  at  a  point,  <{>_,  is  expressed  as  the  sum  of 
its  upward  component,  <P     ,  and  its  downward  component,  A,  .   The  peak 
fraction  corresponding  to  a  given  component  of  bubble  flow  resulting 
from  a  bimodal  cross-covariance  function  is  equal  to  the  corresponding 
time  fraction  of  bubble  flow. 


Experimental 

The  relationship  developed  in  the  preceding  section  between  the 
cross-covariance  peak  fraction  and  the  time  fraction  of  bubble  flow  was 
used  to  analyze  the  bimodal  cross-covariance  functions  characteristic 
of  bubble  flow  in  the  downflow  section  of  the  airlift  tower.   This 
section  describes  the  facilities  and  procedures  used  to  generate  these 
functions. 
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Facilities 

The  description  of  the  experimental  equipment  consists  of  three 
parts:   the  probe,  the  signal  analysis  system  and  the  airlift  tower. 

Probe:   The  presence  of  a  gas  or  liquid  phase  at  two  points  in 
the  flow  system  was  detected  by  a  two  needle  resistivity  probe,  which 
used  the  change  in  electrical  conductivity  between  the  phases  to  detect 
bubbles.   The  probe  design  is  shown  in  Fig.  5.   The  two  needles  were 
fashioned  from  Platinum  10%  Rhodium  wire  with  0.02  cm  O.D.   The  wire 
was  honed  to  a  0.0012  cm  tip  radius  and  inserted  into  a  stainless  steel 
tube  body  with  0.635  cm  O.D.   The  needles  were  insulated  with  a  fast 
drying  varnish  such  that  only  the  tip  was  exposed.   The  tube  served  as 
a  common  return  electrode  for  the  needles. 

The  tube  was  flattened  near  the  needle  end  to  accommodate  the 
needles  and  to  present  the  least  resistance  to  flow.   Fractional 
resistance  to  bubble  flow  was  minimized  by  angling  the  needles  away 
from  the  detection  path  as  shown  in  Fig.  5.   In  the  airlift  tower 
bubbles  were  detected  flowing  in  both  the  upward  and  downward  directions. 
The  needles  were  angled  against  the  downward  direction  of  flow  to  ensure 
the  piercing  and  detection  of  bubbles  moving  downward.   Independent  tests 
in  a  bubble  column  indicated  that,  given  the  same  bubble  flow,  the  probe 
with  needles  angled  against  the  flow,  as  shown  in  Fig.  5,  detected 
25  percent  more  bubbles  than  the  probe  with  needles  angled  away  from  the 

flow. 

(13) 

Herringe  and  Davis    '  used  stainless  steel  suture  needles  for 

resistivity  probes.   Stainless  steel  needles  were  used  in  the  initial 
phases  of  this  work  but  severe  corrosion  problems  made  their  use 
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unacceptable.   Plantinum  10%  Rhodium  alloy  and  Platinum  needles  both  performed 
well,  and  Platinum  10%  Rhodium  was  selected  because  of  its  durability. 

The  distance  between  the  two  needles  of  the  probe  was  determined  by 
photographing  the  probe  tips  with  a  known  reference  length.   A  35  mm  camera 
with  a  50  mm  lens  was  used  with  a  macro  +10  lens  added  to  increase  the  size 
of  the  image  on  the  negative.   Comparative  lengths  were  measured  from  11 
cm  x  17  cm  prints.   The  distance  between  probe  tips  varied  from  probe  to 
probe,  and  ranged  from  0.41  cm  to  0.51  cm. 

Signal  Analysis  System:  The  block  diagram  of  the  resistivity  probe 
system  is  shown  in  Fig.  6.  A  10  kHz  alternating  current  was  supplied  to  the 
common  electrode  by  an  Interstate  Electronics  Corporation  F-74  signal  generator. 
The  conductivity  of  the  X  and  Y  electrical  circuits  changed  instantaneously 
when  bubbles  enveloped  the  respective  probe  tips.   The  resulting  signal  was 
a  10  kHz  current  with  reduced  amplitude  coinciding  with  the  presence  of  a  bubble 
at  the  electrode.   (see  "Raw  Signal"  in  Fig.  5).   This  signal  was  passed 
through  a  comparator  which  grounded  the  signal  with  amplitudes  greater 
than  a  variable  reference  voltage  (see  "After  Comparator"  in  Fig.  5).   The 
frequency  of  this  signal  divided  by  the  input  signal  frequency  yields  the 
fraction  of  time  the  needle  was  enveloped  by  bubbles.   A  Fluke  1900A  Multi- 
counter  was  used  to  obtain  frequencies  of  the  processed  signal  averaged 
over  one  second.   The  average  of  this  value  over  a  40  second  period  of  time 
was  used  to  obtain  the  time  averaged  void  fraction. 

After  the  comparator,  the  signal  was  treated  by  a  frequency  to  voltage 
converter.   The  signals  generated  here  were  square  wave  forms  (see  "After 
Frequency  to  Voltage  Converter"  in  Fig.  5).   These  were  treated  by  a 
Honeywell  Model  SAI-43A  Correlator  and  Probability  Analyzer.   The  analyzer 
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discretized  the  signals  with  a  selected  sample  time  increment.  The  analysis 
was  conducted  in  real  time,  and  over  450,000  observations  were  taken.   The 
sampling  time  was  over  90  seconds  which  was  sufficient  to  ensure  reproduci- 
bility of  the  results. 

Autocovariance  and  cross-covarlance  functions  were  viewed  on  an 
oscilliscope,  and  a  marker  pulse  was  employed  to  locate  the  value  of  h 
corresponding  to  the  covariance  peaks.   The  values  of  h  ranged  from  -200 
to  200,  and  the  sample  time  increment,  At,  was  chosen  in  such  a  way  that  the 
peak  value  of  h  was  near  plus  or  minus  100.   The  analysis  was  conducted  in  the 
clip  mode  of  the  correlator  which  transforms  the  incoming  signals  into  a  +  1 
square  wave.  Since  the  input  to  the  correlator  was  of  the  square  wave  form,  the 
use  of  the  clip  mode  did  not  affect  the  signal  form  and  made  the  most 
efficient  use  of  the  correlator.  Cross  correlation  with  and  without  the 
clip  mode  yielded  identical  time  lag  information. 

(21) 
Airlift  Tower:   The  split  cylinder  airlift  tower     studied  held  a  liquid 

volume  of  10.25  £  and  had  a  15.24  cm  I.D.   The  baffle  height  was  50.  8  cm 

and  the  distance  from  the  lower  edge  of  the  baffle  to  the  bottom  of  the 

tower  was  5.5  cm.   The  ungassed  liquid  was  level  with  the  top  of  the 

baffle.   Air  was  filtered  and  the  flowrate  was  monitered  by  a  Fisher  and 

Porter  (T2-1308/2)  rotameter  which  had  been  calibrated  with  a  wet  test 


Procedures 

Since  the  cross-covariance  analysis  technique  presented  in  this  work  is 
relatively  unknown,  an  investigation  of  the  validity  of  the  technique  for 
slug  bubble  flow  through  a  tube  was  conducted.   The  method  used  to  verify  the 


technique  and  the  procedures  used  in  studying  the  downflow  section  of  the 
airlifts  tower  are  presented. 

Verification  of  Cross-Covariance  Technique:   Slug  bubble  flow  in  a 
0.5  cm  I.D.  latex  tube  was  examined  to  identify  the  cross-covariance 
functions  for  a  simple  unidirectional  flow  system.   A  constant 
liquid  head  and  gas  entrainraent  as  shown  in  Fig.  7  provided  constant 
gas-liquid  flow  conditions.  The  probe  was  inserted  into  a  latex  tubing  and 
the  response  signals  were  recorded  on  a  Hewlett  Packard  tape  recorder. 

The  cross-covariance  function  of  the  response  signals  recorded  on  tape 
was  obtained  for  a  given  period  of  time.   The  tape  direction  was  then 
reversed  and  the  correlation  resumed.   The  resulting  cross-covariance  function 
was  bimodal  and  approximated  that  found  for  two-directional  bubble  flow. 
Since  the  signals  correlated  in  the  forward  and  reverse  directions  were 
identical,  the  time  fraction  of  correlation  conducted  in  the  forward  tape 
direction  corresponded  directly  with  the  time  fraction  of  downward  bubble 
flow.   If  equation  (38)  is  valid,  the  peak  fraction  of  downward  bubble  flow 
will  be  equal  to  the  corresponding  time  fraction  of  flow.   The  autocorrelation 
functions  of  the  response  signals  were  obtained  to  verify  the  random  nature 
of  the  appearance  of  bubbles  at  points  X  and  Y. 

Analysis  of  Bubble  Flow  in  the  Airlift  Tower:   The  airlift  tower  was 
studied  with  an  air-water  system  at  superficial  gas  velocities  ranging 
from  0.0683  m/s  to  0.3115  m/s.   The  liquid  temperature  was  17.4  +  2.7°C 
and  the  air  temperature  was  21°C. 

Open  end  manometers  were  used  to  obtain  pressure  drop  and  regional  hold- 
up data.   These  manometers  were  located  at  the  bottom  and  top  of  each  baffle 
on  the  upflow  and  downflow  sections  of  the  tower.   The  overall  gassed  liquid 
level  was  measured  to  determine  the  overall  gas  hold-up. 
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The  average  liquid  circulation  rate  was  determined  by  hot  water  tracer 

(21) 
analysis     in  the  downflow  section  of  the  tower.   The  responses  from 

Copper-Constantan  thermocouples  centered  in  the  tower  and  43.9  cm  apart  were 

recorded  on  a  Honeywell  Electronite  19  two-pen  stripchart  recorder.   The 

distance  between  peak  centers  was  converted  into  real  time.  The  resulting 

liquid  velocity  is  a  cross-sectional  averaged  quantity  found  from  the 

average  time  of  fluid  travel  between  two  points  in  the  tower.   The  velocity 

changes  inversely  with  the  local  void  fraction.  Therefore,  the  value  of 

interest  is  the  liquid  circulation  rate,  found  as  the  product  of  the  bulk 

liquid  velocity  and  the  average  cross-sectional  area  occupied  by  the  liquid: 

QL  =  vL(l  -  Hd)A  (41) 

2 
where  A  is  the  cross  sectional  area,  m  ,  and  H,  is  the  hold-up  in  the  down- 
flow  section  of  the  airlift  tower. 

The  bubble  flow  was  examined  at  two  levels  in  the  downflow  section  of 
the  airlift  tower.  Level  1  corresponded  to  the  plane  perpendicular  to  the 
baffle  and  4.5  cm  below  the  upper  edge  of  the  baffle.  Level  2  was  located 
1.8  cm  above  the  lower  edge  of  the  baffle.  Both  levels  were  assumed  to  be 
normal  to  the  net  bubble  flow.  At  each  level  the  flow  field  was  examined 
at  four  radial  positions  for  each  of  four  angular  positions.  The  flow  field 
was  examined,  therefore,  at  a  total  of  sixteen  positions  for  each  level. 

Results  and  Discussion 

This  section  contains  two  parts.   The  first  part  examines  the  validity 
of  the  cross-covariance  analysis  technique  described  in  the  theoretical  section. 
The  second  part  presents  and  analyzes  the  results  obtained  by  applying  this 
technique  to  the  airlift  tower. 
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Verification  of  Cross-Covariance  Technique. 

Inherent  in  the  derivation  of  equation  (38)  was  the  assumption  that  the 
presence  of  a  bubble  at  a  point  in  a  flow  field  is  a  random  occurrence.   The 
presence  of  cyclic  phenomena  in  a  time  series  can  be  detected  through  the 
autocovariance  function.   The  autocovariance  function  has  a  maximum  value  at 
a  lag  of  zero  and  is  symmetric  about  zero.   The  presence  of  harmonic  peaks  at 
h  ^  0  indicates  the  presence  of  a  cyclic  process.   The  autocovariance  functions 
for  the  two  electrode  response  signals,  shown  in  Fig.   §,  indicate  that  the 
bubbles  occur  randomly  with  respect  to  time.   The  spread  of  the  peak  at  zero 
lag  is  due  to  the  fact  that  more  than  one  sample  observation  were  taken  for 
each  pulse,  as  discussed  previously.   The  autocorrelation  function  is 
approximated  by  equation  (13a)  with  d  =  0.   A  comparison  of  equation  (13a) 
with  an  experimental  autocorrelation  function  is  presented  in  Fig.  9  with 
£  *  18,   The  agreement  is  good. 

The  cross-correlation  function  is  approximated  by  equation  (13b).   This 
equation  with  d  =  72  and  I  =   20  is  compared  to  an  experimental  cross-correlation 
function  in  Fig.  10.   The  spread  of  the  experimental  cross-correlation  function 
is  greater  than  that  predicted  by  equation  (13b) .   This  can  be  explained  by 
observing  that  the  theoretical  lag,  d,  has  a  constant  value  of  zero  in  the 
case  of  autocorrelation.   This  is  not  true  for  cross-correlation,  and  indeed, 
the  theoretical  lag  for  many  systems  can  be  expected  to  be  a  function  of  time 
with  constant  mean,  as  explained  below. 

The  cross-correlation  function,  by  analogy  with  Eq.  (13),  can  be 
approximated  by 

k  =  1 
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where  d  as  a  function  of  k  is  to  be  determined.   If  d(k)  has  a  skewed 
distribution  the  modal  value  of  p   (h)  will  be  at  a  lag,  h,  between  the 
mean  of  d(k)  and  the  modal  value  of  d(k).   The  time  lag  indicated  by  an 
experimentally  obtained  cross-correlation  function  can  be  considered  as  a 
mean  value  contingent  upon  the  assumption  that  the  mean  time  lag  and  modal 
time  lag  are  close.  The  distribution  of  instantaneous  velocities,  therefore, 
is  assumed  normal.   Since  the  theoretical  lag  as  a  function  of  time  is 
unknown,  the  cross-correlation  function  was  approximated  by  equation  (13b), 
under  the  assumption   that  the  velocity  distribution  is  symmetric  about  the  mean. 
The  theoretical  time  lag  can  be  expressed  as 

d  -  i  (42) 

where   is  the  theoretical  time  lag  in  units  of  time.   The  number  of  obser- 
vations per  pulse  can  be  expressed  similarly  as 


I  -  £  (43) 


where  m   is  the  length  of  time  a  bubble  is  detected  at  a  point.   Equation  (32) 
can  now  be  expressed  as 
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where  P  represents  the  upflow  theoretical  time  lag,  the  downflow  theoretical 
time  lag  or  the  difference  between  them  [see  equation  (32]. 

In  order  that  I   be  greater  than  one  the  constraint  is  added  that 
0  <  m  <_  at.   R  has  a  minimum  value  at  At  =  m  and  a  maximum  value  in  the  limit 
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as  At  approaches  zero.  From  L'Hopitals  rule, 
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(45) 


Since  velocity  can  be  expressed  as 


,  - 1  -  a  (46) 

P   m 


where  B  is  the  distance  between  points  X  and  Y  and  D  is  the  chord  of  the 
bubble  seen  by  the  probe,  the  ratio  P/m  can  be  expressed  as  B/D. 

Four  operational  constraints  now  become  apparent.   The  sample  increment, 
At,  should  be  chosen  large  enough  to  minimize  R,  but  small  enough  to 
maintain  the  sensitivity  of  the  technique  to  signal  variations.  The 
distance  between  points  X  and  Y  should  be  sufficiently  large  to  maximize 
the  ratio  B/D,  but  sufficiently  small  to  ensure  the  correlation  between 
X(t)  and  Y(t). 

Representative  parameter  values  for  slug  bubble  flow  in  a  tube  were 
m  =  1.5  is,  At  ■  0.1  ms  and  P  =  7.9  ms,  and  the  resulting  value  of  R  was 
0.0043.   This  apparently  can  be  neglected  when  compared  to  unity.   The 
assumptions  expressed  as  equation  (32)  are,  therefore,  empirically  justified. 

A  typical  bimodal  cross-covariance  function  found  from  the  slug  bubble 
flow  in  a  tube  (see  Procedures  -  Verification  of  Techniques)  is  presented 
in  Fig.  11.   The  time  fraction  of  upward  bubble  flow  is  0.625.   If  equation 
(38)  is  valid,  the  cross-covariance  peak,  fraction  of  downward  bubble  flow  will 
be  equal  to  the  corresponding  time  fraction  of  flow.  Fig.  12  shows  that 
the  peak  height  fraction  yields  a  reasonable  approximation  of  the  time  fraction 
of  bubble  flow  in  the  corresponding  direction. 
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The  lags  obtained  from  the  bimodal  cross-covariance  function  were 
identical  to  those  obtained  from  the  individual  cross-covariance  functions 
for  unidirectional  flow. 

Analysis  of  Bubble  Flow  in  the  Airlift  Tower 

The  experimentally  determined  gas  hold-ups  in  the  upflow,  downflow  and 
head  regions  of  the  airlift  tower  for  the  air-water  system  are  presented  as 
functions  of  superficial  gas  velocity  based  on  the  upflow  cross-sectional  area 
in  Fig.  13.   The  high  hold-up  in  the  head  region  was  probably  due  to  entrainment 
of  surface  air.   The  overall  hold-up  is  shown  in  Fig.  14  as  a  function  of 
the  superficial  gas  velocity. 

In  Fig.  15  the  cross-sectional  average  linear  liquid  velocity  in  the 
downflow  section  is  shown  as  a  function  of  the  superficial  gas  velocity. 
This  was  determined  from  the  average  time  of  fluid  travel  between  two  points 
in  the  tower  as  described  in  the  Section:   Experimental.   In  Fig.  16  the 
liquid  circulation  rate  is  shown  as  a  function  of  the  superficial  gas  velocity. 
It  can  be  seen  that  the  liquid  circulation  rate  increased  approximately 
linearly  with  the  superficial  gas  velocity. 

The  autocovariance  of  typical  pair  of  response  signals  from  the  downflow 
section  of  the  airlift  tower,  shown  in  Fig.  17,  indicates  that  the  bubble  flow 
in  the  airlift  tower  is  essentially  random.   A  typical  cross-covariance  function 
generated  from  the  downflow  section  of  the  airlift  tower  is  shown  in  Fig.  18. 
The  bimodal  characteristic  of  the  function  is  apparent,  with  one  dominant 
peak  appearing  on  each  side  of  the  zero  lag.   The  lowest  value  of  the  cross- 
covariance  function  occurs  at  a  zero  lag.   This  observation  is  consistent 
with  the  assignment  of  a  zero  value  to  the  cross-covariance  function  at  a 
zero  lag,  as  explained  in  the  theoretical  section. 

The  total  time  averaged  void  fraction  at  a  point  was  found  using  the 
procedure  described  in  the  experimental  section.   The  validity  of  using 
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cross-covariance  peak  fractions  to  estimate  the  time  fraction  of  bubble  flow 
in  a  given  direction  was  shown  in  the  preceding  section.  The  time  fraction,  found 
in  this  manner,  was  used  to  separate  the  void  fraction  into  its  upward  and  down- 
ward components  as 

*T  "  *uP  +  *dn  (48) 

From  the  location  of  the  cross-covariance  peaks  the  positive  and  negative  mean 
lags  were  obtained.   This  results  in  downward  and  upward  velocities.   The 
overall  time-averaged  bubble  velocity,  (vif)  ,  can  be  found  from 

(v<(0  -  v  I       -  v,  tj>  (49) 

o    up  up    dn  dn 

where  v   and  v,   are  positive  quantities  with  direction  indicated  by  the 
subscript.   The  overall  time-averaged  bubble  velocity  at  level  1  is  shown 
as  a  function  of  the  dimensionless  radius  in  Fig.  19.   The  parameters  are  the 
superficial  gas  velocity,  v  ,  and  the  angular  position,  0.  The  bubble 
velocity  and  void  fraction  were  assumed  to  be  zero  at  the  walls.   The  outer 
wall  was  located  at  a  dimensionless  radius  of  one,  and  the  inner  wall  was 
located  at  a  dimensionless  radius  of  0.044  or  0.071,  depending  on  angular 
position.   The  dimensionless  radius  at  the  baffle  wall  was  not  zero  due  to 
the  thickness  of  the  baffle  plate.  The  overall  time-averaged  gas  velocity 
was  downward  at  the  outer  wall  and  upward  near  the  baffle.   The  downward 
liquid  velocity  at  level  1,  therefore,  must  be  greatest  near  the  outer  wall. 
The  gas  phase  flow  field  was  essentially  symmetric  with  respect  to  angular 
position,  normal  to  the  baffle  wall,  9. 

The  void  fraction  at  level  1  is  shown  as  a  function  of  the  dimension- 
less radius  in  Fig.  20.   Again,  the  superficial  gas  velocity  and  angular  position 
are  parameters.   The  largest  void  fraction  is  observed  at  a  dimensionless 
radius  greater  than  0.5.   This  is  due,  in  part,  to  the  inability  of  the  probe 
to  detect  all  rising  bubbles  as  explained  previously.   The  separation  of  the 
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total  void  fraction  into  its  upward  and  downward  components  should  be  viewed, 
therefore,  as  leading  to  a  correction  of  downward  bubble  flow  rate  measurements. 
The  upward  bubble  flow  rate  can  be  determined  through  mass  balance  over  the 
downflow  region. 

The  void  fraction,  $_,  can  be  found  as 

*i  "!</«/  vdrde  (50) 

where  integration  is  performed  over  the  plane  normal  to  the  flow.  The 
downward  and  upward  flow  rates  of  air  through  that  plane  are  found, 
respectively,  as 

Q.   =  f     /R  v,  4  ,   r  dr  d6  (51a) 

^A,    o   o    dn  dn 
dn 

Q,    =  f     /R  v  4    r  dr  d9  (51b) 

A     o   o    up  up 

up 

Integration  of  equations  (50)  and  (51)  was  performed  graphically  through  the 

use  of  a  planimeter.  Values  of  v,  4,  r,  v  4  r  and  4  r  are  presented 
r  dnTdn    up  up       T 

as  functions  of  the  dimensionless  radius  with  the  angular  position,  6,  as 

a  parameter,  in  Figs.  21,  22  and  23,  respectively. 

The  integral  of  v.  A,  r  with  respect  to  r  is  shown  as  a  function  of 
■       dn  dn 

the  angular  position  in  Fig.  24.   Maxima  are  observed  at  an  angle  of  0.5tt 

corresponding  to  the  line  perpendicular  to  the  baffle  wall.  In  Fig.  25 

the  integral  of  the  upward  velocity,  v  4   r,  with  respect  to  r  is  presented 

as  a  function  of  the  angular  position.   This  function  appears  to  be  bimodal 

with  maximum  values  at  9  ■  0.2tt  and  0.8tt  and  local  minima  at  9  =  0,  0.5tt 

and  it.  The  integral  of  4  r  with  respect  to  r  is  shown  as  a  function  of 

the  angular  position  in  Fig.  26.   The  values  presented  Fig.  24,  25  and  26  were 
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graphically  integrated  to  obtain  the  downward  air  flow  rate,  upward  air  flow 
rate  and  void  fraction  at  level  1.   A  similar  analysis  was  conducted  at 
level  2,  located  1.8  cm  above  the  lower  edge  of  the  baffle.   The  conducti- 
vity probe  was  able  to  detect  bubble  flow  only  at  the  highest  superficial 

gas  velocity  studied,  0.3115  m/s.   Values  of  v  <t  r,  v.  A,  r  and  A_r  are 
&         *  *  up  up    dnrdn       T 

presented,  respectively,  in  Figs.  27,  28  and  29.   The  integrals  of  these 
quantities  with  respect  to  the  radius  are  shown  as  functions  of  the  angular 
position  in  Figs.  30  and  31. 

The  gas  flow  rates  so  determined  are  presented  in  Fig.  32.   The  downward 
flow  of  air  at  level  1  with  a  superfical  gas  velocity  of  0.3115  m/s  based 
on  the  upflow  cross-section  area  was  2.88  x  10   m  /s .   This  is  10.2  percent 
of  the  incoming  air  flow  rate.   The  downward  flow  of  air  at  level  2  can  be 

considered  as  the  flow  of  air  carried  into  the  upflow  section  by  the  liquid 

-4   3 
circulation.   This  was  found  to  be  0.2508  x  10   m/s  (0.9  percent  of  the 

incoming  gas  flow  rate).   A  mass  balance  on  the  downflow  section  at  the 

superficial  gas  velocity  of  0.3115  m/s  yielded  the  gas  flow  rate 

in  the  upward  direction  at  level  1  as  2.63  x  10   m/s.   The  value  determined 

through  integration  was  0.875  x  10   m/s,  and  therefore,  only  one  third 

of  the  bubbles  moving  in  the  upward  direction  was  detected. 

It  is  important  to  note  that  uncertainty  is  present  in  the  integrated 

results  due  to  the  paucity  of  data  points  obtained  at  each  level.   The  void 

fractions  found  through  integration  at  levels  1  and  2  are  compared  in 

Fig.  33  to  the  gas  hold-up  determined  through  pressure  drop  measurement 

in  the  downflow  region.   In  this  integration  the  void  fraction  was  assumed 

to  be  zero  at  the  wall  as  previously  discussed.   If  a  non-zero  void  fraction 

at  the  walls  is  assumed  the  average  void  fraction  becomes  larger  by  30  to 

50  percent.   These  results  are  also  presented  in  Fig.  33  as  functions  of 

superficial  gas  velocity.   It  can  be  noted  that  the  void  fraction  at 
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level  1  is  much  larger  relative  to  the  gas  hold-up  at  low  superficial 
gas  velocities  than  it  is  at  higher  superficial  gas  velocities.   This 
is  expected,  since  the  region  in  which  a  larger  concentration  of  gas 
bubbles  is  observed  increases  with  superficial  gas  velocity.   The  void 
fractions  presented  in  this  work  are  low  since  two  thirds  of  the  bubbles 
moving  in  upward  flow  were  not  detected  due  to  the  bias  of  the  probe  to 
detecting  bubbles  in  the  downward  flow. 
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Conclusions 

Bubble  flow  in  the  downflow  section  of  the  airlift  tower  has  net  upward 
and  downward  components.   A  technique  based  on  time  series  analysis  was 
developed  to  obtain  the  velocity  and  time  fraction  of  bubble  flow  for  each 
component.   Resistivity  probes  were  employed  to  obtain  dependent  time  series 
indicating  the  presence  of  the  gas  or  liquid  phase  at  a  point.   The  cross- 
covariance  function  of  these  series  was  bimodal  with  local  maxima  giving 
rise  to  the  time  lag  characteristics  of  each  component  of  flow.  The  time 
fraction  of  flow  was  found  to  be  equal  to  the  peak  fraction  of  the  bimodal 
cross-covariance  function,  linearly  translated  so  that  the  function  at  zero 
lag  has  a  value  of  zero.   This  technique  was  developed  theoretically  and 
verified  experimentally.   A  number  of  significant  conclusions  can  be 
drawn  concerning  the  technique. 

1)  The  appearance  of  bubbles  at  a  point  in  a  gas-liquid  flow  field  is 
a  random  occurence. 

2)  The  autocorrelation  function  for  random  bubble  flow  is  modeled  by 
equation  (13) . 

3)  The  cross-correlation  function  is  distinguished  from  the  auto- 
correlation function  in  that  the  theoretical  lag  may  be  non  zero  and  is  not 
constant.  A  reasonable  model  for  the  cross-correlation  function  is 
presented  as  equation  (41) . 

The  technique  developed  was  employed  in  a  study  of  the  downflow  section 
of  the  airlift  tower.   The  significant  conclusions  are: 

1)  Bubble  flow  in  the  downflow  section  of  the  airlift  tower  has  upward 
and  downward  components  of  velocity. 

2)  The  downflow  of  air  at  a  level  of  4.5  cm  below  the  top  edge  of  the 

-4   3 
baffle  ranged  from  1.23  x  10   m  /s  at  a  superficial  gas  velocity  of  0.0683  m/s 

to  2.88  x  10   m  /s  at  a  superficial  gas  velocity  of  0.3115  m/s. 
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3)  The  flowrate  of  air  carried  from  the  downflow  section  into  the 

-4  3 
upflow  section  was  0.2508  x  10   m  /s  at  a  superficial  gas  velocity  of 

0.3115  m/s. 

4)  The  bubble  flow  field  was  essentially  symmetric  about  the  angular 
position. 

5)  The  bubble  flow  at  the  upper  level  of  the  airlift  tower  was  downward 
at  a  large  radius  and  upward  near  the  center  of  the  tower. 

The  velocities  and  time  fractions  of  flow  in  complex  two  phase  flow 
systems  can  be  found  through  an  analysis  of  the  cross-covariance  functions 
characteristic  of  the  flow  field.   The  techniques  developed  in  this  work 
to  extract  this  information  from  bimodal  cross-covariance  functions  are 
general,  and  should  be  employed  in  other  multi-phase  systems. 
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Notation 

A 

Cross-Sectional  area 

B 

Distance  between  points  X  and  Y 

D 

Chord  of  bubble  seen  by  probe 

d 

Theoretical  lag 

E 

Expectation 

Hd 

Hold-up  in  the  downflow  section  of  the  airlift  tower 

h 

Lag 

k 

Sample  number 

K 

Total  number  of  observations 

L 

Height  of  cross-covariance  function  peak  relative  to  the  value  of 

cross-covariance  function  at  zero  lag 

I 

Number  of  observations  per  pulse 

m 

Length  of  time  the  gas  phase  occupies  a  point,  s 

0 
P 

Order  of  probability 

P 

Theoretical  time  lag,  s 

*A 

3 
Flow  rate  of  gas  phase,  m  /s 

^L 

3 
Flow  rate  of  liquid  phase,  m  /s 

r 

Radial  position 

t 

Time,  s 

At 

Sample  increment ,  s 

V 

Local  bubble  velocity,  m/s 

VL 

Bulk  liquid  velocity,  m/s 

V 

s 

Superficial  gas  velocity,  m/s 

W(t) 

Stationary  time  series 

X 

Upper  probe  point  position 

X(t) 

Time  series  corresponding  to  void  fraction  at  X 

Y 

Lower  probe  point  position 
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Y(t)  Time  series  corresponding  to  void  fraction  at  Y 

Z(t)  Stationary  time  series 

Greek 

YyyCh)  Cross  covariance  function  of  X(t)  and  Y(t) 

Y  (h)  Autocovariance  function  of  X(t) 
XX 

8  Angular  position,  rad. 

u  Mean 

pyyCh)  Cross-correlation  function  of  X(t)  and  Y(t) 

p^yCh)  Estimated  cross-correlation  function  of  X(t)  and  Y(t) 

Py^Ch)  Autocorrelation  function  of  X(t) 

p   (h)  Estimated  autocorrelation  function  of  X(t) 

T  Time  lag,  s 

<(>  Local  time  averaged  void  fraction 

4>  Integrated  void  fraction  for  a  given  level 

X  Time  fraction  of  flow 

CI  Peak  fraction 

Subscripts 

dn  Downward  component 

T  Total 

up  Upward  component 

W  Corresponding  to  W(t) 

X  Corresponding  to  X(t) 

Y  Corresponding  to  Y(t) 
Z  Corresponding  to  Z(t) 
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b)   Autocorrelation  Function  for  W(k) . 
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APPENDIX  A:   Construction  of  Conductivity  Probes 
The  presence  of  a  gas  or  liquid  phase  at  a  point  can  be  detected  by 
observing  the  difference  in  conductivity  between  the  two  phases.   The 
probe  used  to  detect  this  difference,  as  shown  in  Fig.  7,  employs  two 
sharp  needles  and  a  stainless  steel  tube  as  electrodes.   The  stainless 
steel  tube  serves  both  as  a  housing  for  the  needle  electrodes,  and  as  a 
common  electrode.   The  needles  are  insulated  except  at  the  tips,  and  a 
10  kHz  alternating  current  is   allowed  to  flow  between  the  common  electrode 
and  the  two  needle  electrodes.   As  shown  in  Fig.  7,  when  a  bubble  envelops 
an  electrode,  the  resistance  to  current  increases  instantaneously,  and  the 
amplitude  of  the  response  signal  from  the  electrode  decreases.   The  time 
interval  for  which  the  signal  amplitude  is  decreased  is  the  time  interval  in 
which  the  bubble  "insulates"  the  electrode. 

Design  Factors 

Four  major  factors  need  to  be  considered  in  the  design  of  a  two-needle 
conductivity  probe. 

1)  The  two  needles  and  the  tube  body  need  to  be  electrically 
insulated.   This  also  means  that  the  probe  must  be  water  tight. 

2)  The  probe  should  be  capable  of  detecting  the  bubbles  which  enter 
the  detection  path.   The  needles  should  be  sharp  enough  to  pierce  the 
bubbles  which  come  into  contact  with  it.   The  best  piercing  action  occurs 
when  the  surface  area  initially  contacting  the  bubble  is  minimized.  A 
suitable  probe  configuration  is  shown  in  Fig.  Ala. 

The  total  surface  area  of  the  needles  contacting  the  bubble  should  be 
minimized  to  reduce  the  effect  of  the  probe  on  bubble  velocity.   The  drag 
exerted  by  the  probe  on  the  bubble  is  greater  for  the  configuration  shown 
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in  Fig.  Ala  than  for  the  configuration  shown  in  Fig.  Alb.  The  probe  shown  in 
Fig.  Alb,  however,  will  not  pierce  the  bubbles  effectively.   A  compromise 
configuration,  shown  in  Fig.  Ale,  is  chosen  to  provide  good  piercing  and  a 
small  amount  of  drag . 

3)  The  needles  should  maintain  their  geometric  integrity.   The  distance 

between  the  probe  tips  should  be  constant,  and  the  material  used  should  be 

(13) 
corrosion  resistant.  Herringe  and  Davis     used  stainless  steel  surgical 

needles  as  electrodes;  however,  severe  corrosion  problems  were  encountered 

when  they  were  used  in  the  present  study.   The  large  current  density  at  the 
needle  tips  makes  the  use  of  a  more  noble  metal  desireable.   Both  Platinum 

and  Platinum  10%  Rhodium  performs  well.  The  platinum  alloy  is  selected 

because  of  its  superior  resistance  to  deformation. 

4)  The  distance  between  the  tips  of  the  needles  should  be  large  enough 
to  be  accurately  measured,  yet  small  enough  that  a  bubble  contacting  one 
needle  could  be  expected  to  contact  the  other.   For  the  probe  described  here 
the  distance  ranges  from  0.41  to  0.51  cm.  The  distance  can  be  measured  by 
photographing  the  probe  tips  along  with  a  micrometer  which  provided  a 
reference  length. 


Procedure 

1)   Preparation  of  tube  body.   In  the  present  probe  the  tubular  body  is 
fashioned  from  a  stainless  steel  tube  with  a  0.635  cm  0.D,  0.508  cm  I.  D. 
length,  approximately  15  cm  long,  is  cut  from  a  straight  section  of  tubing. 
The  length  of  the  tube  used  is  dependent  upon  the  size  of  the  system  for 
which  the  probe  is  intended. 

As  shown  in  Fig.  A2a,   an  electrical  connection  is  attached  to  the  tube 
at  one  end.   This  is  done  by  wrapping  the  wire  around  the  tube  and  soldering 
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it,  or  by  drilling  a  hole  in  the  tube  and  inserting  the  wire  into  it  before 
soldering.   A  resin-core  solder  with  a  zinc  chloride  flux  is  used  for  all 
connections. 

The  compression  fitting  is  placed  on  the  tube  with  the  threaded  end 
facing  the  end  of  the  tube  that  is  not  used  for  the  electrical  connection. 
This  end  was  placed  in  a  vise  and  the  last  centimeter  of  the  tube  was 
pinched  slightly.   The  pinching  of  the  tube,  shown  in  Fig.  A2a  gives  the 
end  of  the  tube  an  elliptical  cross-section.   This  shape  allowed  the  needle 
to  be  placed  farther  apart  than  they  could  be  with  an  unflattened  tube.   For 
a  given  distance  between  needles  a  smaller  diameter  tube  can  be  used  for  the 
body  if  the  end  is  flattened.   This,  along  with  the  streamlining  effect  of 
the  flattened  section,  minimizes  the  disturbance  of  flow. 

2)  Construction  of  needle  support  assembly.   The  needle  support  tubes 
are  assembled  as  a  separate  unit  which  is  then  inserted  into  the  tube  body. 
The  needle  support  tubes,  with  0.165  cm  O.D.  and  0.089  cm  I.  D.,  are  2.5  cm 
in  length.   A  33  degree  bend  is  formed  0.5  cm  from  the  end  of  each  support 
tube.   The  exact  angle  is  not  critical  as  long  as  both  tubes  have  the  same 
angle. 

The  electrical  leads  are  attached  to  the  straight  ends  of  these  tubes. 
This  is  done  by  coating  the  wires  with  molten  solder  and  inserting  them  into 
the  support  tubes.  The  two  support  tubes  are  clamped  by  the  connected  wires 
in  a  table  vise.   Quick  setting  epoxy  is  used  to  join  the  support  tubes 
together  and  insulate   them  from  each  other.   The  finished  product,  shown  in 
Fig.  A2b,  should  be  small  enough  to  fit  into  the  tube  body. 

3)  Joining  of  needle  support  tube  assembly  and  tube  body.   The  electrical 
leads  to  the  support  tube  assembly  were  threaded  into  the  tube  body,  and  the 
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assembly  is  Inserted  into  the  elliptical  end  of  the  tube  body  as  shown  in 
Fig.   A2c.  Epoxy  is  used  to  join  the  support  tube  assembly  to  the  tube  body 
to  avoid  leakage  of  water  into  the  probe,  and  to  insulate  the  support  tubes 
from  the  tube  body. 

4)  Insertion  of  Platinum  10%  Rhodium  needles.   Platinum  10%  Rhodium 
0.051  cm  0.  D.  wire  is  honed  to  a  0.0012  cm  tip  radius  to  form  the  needle 
electrodes.   The  needle  support  tubes  are  filled  with  molten  solder  and 
the  platinum  alloy  needles  inserted.  A  thin  layer  of  quick-setting  epoxy 
is  applied  to  the  exposed  needle  support  tubes  and  to  the  solder  connection 
The  result  is  shown  in  Fig.  A2d.   The  use  of  support  tubes  in  the  probe 
assembly  makes  the  exchange  of  new  needles  for  damaged  ones  an  easy  process. 

5)  Insulation  of  needle  electrodes.  The  final  step  in  the  construction 
of  the  conductivity  probe  is  to  spray  the  probe  needles  with  Krylan  Crystal 
Clear  Protective  Spray  Coating  No.  1302.   This  coating  insulates  the  needles. 
The  surface  tension  is  large  enough  that  the  sharp  tips  of  the  needles 
remained  uncoated.   As  needed,  the  spray  coating  is  reapplied  during  use. 
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Fig.  Al.   The  Effectiveness  of  Probe  Configurations; 

a)  Good  Piercing  Action,  Large  Amount  of  Drag, 

b)  Poor  Piercing  Action,  Small  Amount  of  Drag, 

c)  A  Compromise  Between  a)  and  b) . 
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Fig.  A2.   Schematic  Diagram  of  Systematic  Conductivity  Probe 
Construction; 

a)  Preparation  of  Tube  Body, 

b)  Construction  of  Needle  Support  Tube  Assembly, 

c)  Joining  of  Support  Tube  Assembly  and  Tube  Body, 

d)  Insertion  of  Platinum  10%  Rhodium  needles. 


A.   Electrical  Connection  to  Tube  Body,  B.   Compression 
Fitting,  C.   Flattened  End  of  Tube  Body,  D.   Copper 
Needle  Support  Tubes,  E.   Epoxy,  F.   Electrical  Connections 
to  Support  Tube,  G.   Platinum  10%  Rhodium  needles. 
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APPENDIX  B.   Block  Diagram  of  Response  Signal  Treatment  System 
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APPENDIX  C.   Data  Summary  for  Airlift  Tower 


Table  CI.   Physical  Description  of  Airlift  Tower 
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Liquid  Volume 

10.25  I 

Diameter 

15.2  cm 

Baffle  Height 

51.0  cm 

Baffle  Width 

0.635  cm 

Distance  from  Bottom  of  Column  to 

Bottom  Edge  of  Baffle 

5.5  cm 

Sparger  Tube  Height 

5.5  cm 

Sparger  Tube  I.D. 

0.95  cm, 

Upflow  Cross-Sectional  Area 

90.7  cm"" 

Level  1  Position 

52.0  cm 

Level  2  Position 

7.3  cm 
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Table  C7. 

Intermediate 

Integrated 

Values 

Level 

Angular 

V 

L  v     4 

rdr 

f  v  ,   $  .   rdr 

R 

L   *,rdr 

Position, 
rad. 

(m/s) 

0     up  up 
(m3/s) 

0     dn  dn 
(m3/s) 

0   TT, 
m 

1 

0.2* 

0.0683 

3.88  x 

io-5 

5.10  x  10"5 

2.79   x  10"4 

1 

0.2rr 

0.1412 

2.93 

4.15 

2.52 

1 

0.2* 

0.2307 

3.43 

5.73 

2.88 

1 

0.2ii 

0.3115 

4.72 

7.50 

3.37 

1 

0.4i7 

0.0683 

2.01 

4.01 

2.70 

1 

0.4ir 

0.1412 

2.37 

3.43 

2.23 

1 

0.4it 

0.2307 

1.45 

9.24 

2.36 

1 

0.4* 

0.3115 

1.09 

14.32 

3.02 

1 

0.6n 

0.0683 

1.63 

7.07 

3.50 

1 

0.6rc 

0.1412 

1.83 

4.75 

3.03 

1 

0.6* 

0.2307 

2.99 

14.33 

4.49 

1 

0.6n 

0.3115 

2.46 

14.26 

3.84 

1 

0.3" 

0.0683 

1.83 

3.21 

1.96 

1 

0.8" 

0.1412 

2.60 

3.88 

2.32 

1 

0.8" 

0.2307 

4.24 

7.35 

3.34 

1 

0.8" 

0.3115 

5.33 

9.60 

3.79 

2 

0.2" 

0.3115 

0.395 

1.960 

1.753 

2 

0.4" 

0.3115 

0.703 

0.657 

1.180 

2 

0.6" 

0.3115 

0.157 

0.521 

0.671 

2 

0.8" 

0.3115 

0.155 

0.791 

0.797 
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Table  C8.  Integrated  Flowrate  and  Void  Fraction 


Level 

V 

s 

QA  C*3 
up 

's> 

Q.  (">3/s) 
A. 
an 

*T 

1 

1 
1 

1 

2 

0.0683 
0.1412 
0.2307 
0.3115 

0.3115 

0.602  x 
0.632 
0.780 
0.875 

0.0898 

10"* 

1.23  x 
1.03 
2.35 
2.88 

0.2508 

ID"4 

0.0773 
0.0742 
0.0895 
0.0974 

0.0306 
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Chapter  V 

CONCLUDING  REMARKS 

This  chapter  summarizes  the  conclusions  presented  in  Chapters  III 
and  IV.   Recommendations  for  future  work  are  presented. 

Conclusions 

The  comparison  of  one  and  two-stage  airlift  towers  based  on  mass 
transfer  rates  and  efficiencies  yielded  the  following  conclusions: 

1)  The  use  of  a  multi-stage  airlift  tower  has  merit  in  gas-liquid 
contacting.   Oxygen  transfer  coefficients  (K.  a)  of  over  25  min 

(1500  hr  )  were  obtained  at  a  superficial  gas  velocity  of  2728  cm/min 
for  a  sodium  sulfite-air  system  in  a  two-stage  tower.   At  a  superficial 
gas  velocity  of  2262  cm/min  the  IC  a  for  the  two-stage  tower  was  54 
percent  higher  than  in  the  one-stage  system. 

2)  The  use  of  a  multi-stage  airlift  tower  requires  more  reactor 
volume  than  the  single  stage  tower  of  equal  liquid  volume,  but  the 
oxygen  transfer  coefficient  based  on  dispersion  volume  is  still  27  percent 
higher  than  in  the  one  stage  tower  at  a  superficial  gas  velocity  of  2262 
cm/min. 

3)  At  low  gas  velocities  the  two-stage  and  one-stage  towers  have 
equivalent  mass  transfer  characteristics.   At  a  superficial  gas  velocity 
of  1300  cm/min  the  JL  a  levels  off  for  the  one-stage  tower,  while  the  IC  a 
for  the  two  stage  tower  increases  linearly  with  superficial  gas  velocity 
in  the  range  studied. 

4)  The  two-stage  tower  becomes  more  efficient  than  the  single  stage 
tower  at  oxygen  transfer  rates  larger  than  180  m  mole/S.  hr. 
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Bubble  flow  in  the  downflow  section  of  the  airlift  tower  has  net  upward 
and  downward  components.   A  technique  based  on  time  series  analysis  was 
developed  to  obtain  the  velocity  and  time  fraction  of  bubble  flow  for  each 
component.   Resistivity  probes  were  employed  to  obtain  dependent  time  series 
indicating  the  presence  of  the  gas  or  liquid  phase  at  a  point.   The  cross- 
covariance  function  of  these  series  was  bimodal  with  local  maxima  giving 
rise  to  the  time  lag  characteristics  of  each  component  of  flow.   The  time 
fraction  of  flow  was  found  to  be  equal  to  the  peak  fraction  of  the  bimodal 
cross-covariance  function,  linearly  translated  so  that  the  function  at  zero 
lag  has  a  value  of  zero.   This  technique  was  developed  theoretically  and 
verified  experimentally.   A  number  of  significant  conclusions  can  be 
drawn  concerning  the  technique. 

1)  The  appearance  of  bubbles  at  a  point  in  a  gas-liquid  flow  field  is 
a  random  occurrence. 

2)  The  autocorrelation  function  for  random  bubble  flow  is  modeled  by 
equation  (13). 

3)  The  cross-correlation  function  is  distinguished  from  the  auto- 
correlation function  in  that  the  theoretical  lag  may  be  non/zero  and  is  not 
constant.   A  reasonable  model  for  the  cross-correlation  function  is 
prescribed  as  equation  (41) . 

The  technique  developed  was  employed  in  a  study  of  the  downflow  section 
of  the  airlift  tower.   The  significant  conclusions  are: 

1)  Bubble  flow  in  the  downflow  section  of  the  airlift  tower  has  upward 
and  downward  components  of  velocity. 

2)  The  downflow  of  air  at  a  level  of  4.5  cm  below  the  top  edge  of  the 

baffle  ranged  from  1.23  x  10   m  /s  at  a  superficial  gas  velocity  of  0.0683  m/s 

-4  3 
to  2.88  x  10   m  fa   at  a  superficial  gas  velocity  of  0.3115  m/s. 
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3)  The  flowrate  of  air  carried  from  the  downflow  section  into  the 

-4  3 
upflow  section  was  0.2508  x  10   m  /s  at  a  superficial  gas  velocity  of 

0.3115  m/s. 

4)  The  bubble  flow  field  was  essentially  symmetric  about  the  angular 
position. 

5).  The  bubble  flow  at  the  upper  level  of  the  airlift  tower  was  downward 
at  a  large  radius  and  upward  near  the  center  of  the  tower. 

The  velocities  and  time  fractions  of  flow  in  complex  two  phase  flow 
systems  can  be  found  through  an  analysis  of  the  cross-covariance  functions 
characteristic  of  the  flow  field.   The  techniques  developed  in  this  work 
to  extract  this  information  from  bimodal  cross-covariance  functions  is 
general,  and  should  be  employed  in  other  multi-phase  system. 

Recommendations 

The  scale-up  of  one-stage  and  multi-stage  airlift  towers  needs  to  be 
examined.   The  magnitude  of  improvement  of  mass  transfer  characteristics 
observed  with  20  4  staged  airlift  towers  may  be  different  than  that 
observed  in  larger  systems. 

The  motivations  for  staging  in  fermentation  systems  are  to  increase 
mass  transfer  rates  and  to  achieve  greater  cell  yields  by  employing 
partial  mixing.   The  one-stage  and  two-stage  towers  should  be  examined 
under  continuous  fermentation.   The  organisms  chosen  for  this  study  should 
have  high  growth  rates  with  large  oxygen  consumption  rates.   The  effect 
of  staging  on  mass  transfer  is  more  apparent  under  large  oxygen  transfer 
rates. 

Some  preliminary  studies  of  the  one-stage  tower  under  batch 
fermentation  were  conducted.   The  tower  dimensions  are  described  in  Chapter 
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III.   A  mixed  culture,  obtained  from  a  soil  sample,  was  grown  on  a  glucose 
and  salts  medium  as  shown  in  Table  1.   The  fastest  growing  organisms  in 
the  soil  sample  were  selected  for  use. 

A  paramagnetic  oxygen  analyzer  was  used  to  measure  the  oxygen 
concentration  of  the  exit  gas.  The  dissolved  oxygen  concentration  was 
measured  in  the  upflow  and  downflow  sections  of  the  airlift  tower.   An 
oxygen  balance  yielded  the  values  of  the  oxygen  transfer  coefficient 
which  is  shown  in  Fig.  1  as  a  function  of  cell  concentration  for  the  one- 
stage  airlift  tower.   The  surface  tension  was  measured  to  be  47.5  dynes/cm, 
and  the  liquid  density  was  1.015  g/cm  .   Foaming  problems  were  controlled 
by  use  of  Dow  Corning  silicone  Antifoam  agent.   The  culture  was  in  the 
exponential  growth  phase.   This  technique  was  workable  in  the  20  8.  system. 

Turbulent  two-phase  flow  can  be  characterized  by  its 

(2) 
statistical  properties.   T.  Y.  Chen    observed  bimodal  cross-covariance 

functions  in  solid  flow  through  a  tube.   The  technique  developed  in 

Chapter  IV  for  the  analysis  of  gas-liquid  flow  can  be  applied  to  gas-solid 

and  liquid-solid  flow. 

The  airlift  tower  should  be  investigated  more  thoroughly,  increasing 
the  number  of  data  points  obtained  at  each  level  and  increasing  the  number 
of  levels  investigated.   The  statistical  methods  developed  can  be  employed 
in  the  investigation  of  larger  airlift  towers. 

The  bubble  flow  was  analyzed  in  this  study  along  the  net  upward  and 
downward  directions  of  flow.   The  turbulent  nature  of  the  flow  is  such 
that  a  horizontal  component  of  flow  exists.   This  can  be  examined  using 
similar  techniques. 
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Table  1.   Composition  of  Fermentation  Medium 

C6H12°6  30   8/4 

Nutrient  Broth  3   g/S. 

Na2HP04  3.8  g/£ 

KH2P04  3.5  g/4 

(NH4)2S04  3.5  g/l 

MgS04-7H,0  1   g/i 

Trace  Metal  Solution  3  ml/I 

Trace  Metal  Solution 


27      g/a 

2        g/i 
2        g/* 

2        g/l 

0.8   g/S. 

100     ml/it 


FeCl  • 

6H  0 

Na2MoC 

l4-2H20 

CoS04' 

7H,0 

ZnSO, • 

4 

7H20 

CaCl2 

H3B03 
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One  of  the  more  promising  contactor  designs  for  high  mass  transfer 
rates  with  low  power  consumption  is  the  airlift  tower.  The  mass  transfer 
coefficients^   reported  in  the  literature  for  the  airlift  tower  were 
compared  to  the  bubble  column,  agitated  vessel,  staged  bubble  column  and 
stage  agitated  vessel.   Staging  has  been  shown  to  improve  the  mass  transfer 
characteristics  of  bubble  columns  and  agitated  vessels. 

The  oxygen  transfer  characteristics  of  a  two  stage  split  cylinder  air- 
lift tower  were  investigated  and  compared  to  those  of  a  similar  single 
stage  airlift  tower  of  equal  liquid  volume  using  a  sodium  sulfite-air 
system.   At  superficial  gas  velocities  from  0.120  to  0.200  m/s  no 
difference  in  K^a  was  apparent.   The  (La  was  significantly  larger  in 
the  two  stage  tower  for  a  gas  velocity  between  0.200  and  0.455  m/s. 
At  0.455  m/s  a  K^a  of  25.2  min   was  achieved  in  the  two  stage  system, 
and  at  0.377  m/s  the  two  stage  tower  had  a  54  percent  larger  K. a  than  the 
single  stage.   A  comparison  of  dispersion  volume  based  K  a  showed  a  27 
percent  larger  value  at  a  gas  velocity  of  0.377  m/s.   The  performance 
ratio  for  the  two-stage  tower  was  larger  than  that  for  the  single  stage 
tower  at  oxygen  transfer  rates  greater  than  180  m  male/Jt  hr. 

Bubble  flow  in  the  airlift  tower  has  both  upward  and  downward 
components.   A  technique  based  on  time  series  analysis  was  developed  to 
obtain  the  velocity  and  time-averaged  void  fraction  for  each  component 
of  bubble  flow.   This  technique  was  derived  theoretically  and  verified 
experimentally . 

The  technique  was  employed  in  a  study  of  the  downflow  section  of  an 
airlift  tower  for  an  air-water  system.   The  downflow  of  air  at  the  upper 
level  of  the  tower  ranged  from  1.23  x  10~4  m3/s  at  a  superficial  gas 
velocity  in  the  upflow  section  of  0.0633  m/s  to  2.88  x  10"4  m3/s  at  0.3115  m/s. 


The  flow  rate  of  air  carried  from  the  downflow  section  into  the  upflow 

-4   3 
section  was  0.2508  x  10   m  /s  at  a  superficial  gas  velocity  of  0.3115  m/s. 


